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Abstract 



Based on a microscopic density functional theory we investigate the mor- 
^ I phology of thin hquidlike wetting films adsorbed on substrates endowed with 

c/3 ' well-defined chemical heterogeneities. As paradigmatic cases we focus on a 



single chemical step and on a single stripe. In view of applications in microflu- 
idics the accuracy of guiding liquids by chemical microchannels is discussed. 
Finally we give a general prescription of how to investigate theoretically the 



^ , wetting properties of substrates with arbitrary chemical structures 
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I. INTRODUCTION 

A large variety of experimental techniques has emerged which allows one to endow solid 
substrates with stable, persistent, and well-defined lateral patterns of geometrical or chemical 
nature or a combination thereof ||l|-§. The spatial extension of these regular structures 
ranges from the /im scale down to nm. A major part of the application potential of such man- 
made surfaces is based on their exposure to fiuids and their ability to imprint permanent, 
well-defined structures on adjacent soft matter encompassing simple anorganic and organic 
liquids as well as complex fiuids such as polymer solutions or colloidal suspensions (see, e.g., 
Ref. [0). 

Here we are interested in planar, chemically structured surfaces. Recent experiments 
have explored the wetting and adsorption properties of fiuids and polymer solutions on 
such heterogeneous surfaces |]8|-[15|; some studies even have demonstrated the feasibility to 
control the growth of biological systems by attaching them to structured surfaces P^JTT]] and 
to recognize biological molecules, e.g., proteins, selectively by bringing them into contact 



with nanostructured surfaces [|l8l , p!9| . In the context of microfiuidics pO| , pT| one is interested 
in guiding small quantities of valuable liquids to designated sites where they can either be 
analyzed or undergo chemical reactions. This transportation problem can be solved by using 
microgrooves or microchannels. Alternatively, liquids can be guided along chemical lanes 
on fiat substrates. Such an integrated network can be used to build chemical chips and 
tiny chemical factories (see, e.g., Ref. [^). This development has been greatly facilitated 
by the emergence of microcontact printing p3|-[26[] as an important technique to cover solid 
substrates with designed chemical patterns. With this technique one can create surface 
patterns which consist of regular monolayer patches composed of different chemical species. 
These patterns are anchored at a homogeneous and flat substrate such that the resulting 
decorated surface remains flat on the molecular scale. Such a patchwork of, e.g., hydrophilic 
and hydrophobic areas with lateral extensions in the sub-micrometer range can be designed 
with high precision and turns out to be rather robust. 

In view of the desired miniaturization of such structures the question arises to which 
extent thermal fluctuations limit the ability to keep adsorbed liquids with high lateral pre- 
cision on the designated chemical patterns without spilling. Moreover, one would like to 
know which chemical and structural architecture of the man-made surface is favorable for 
miniaturization. This requires to decode the relation between the structural properties of the 
highly inhomogeneous liquids and the molecular interaction potentials of the fluid particles 
and of the various, speciflcly arranged substrate particles. 

We address some of these issues for the case of thermal equilibrium by suitable tools 
of statistical mechanics. Different theoretical models have already been applied to study 
various basic properties of adsorption of liquids on structured substrates. Exact virial the- 
orems and compressibility sum rules, as they can be formulated for fluids in contact with 



homogeneous walls, have been derived also for heterogeneous substrates |27|. In this context 



phenomenological interface displacement and square-gradient models ||Tl|,pB|-Bl] as well as 



lattice gas models p5|-p7[] highlight the general behavior of liquidlike wetting layers and their 



singular properties at surface or morphological phase transitions. More sophisticated the- 
ories such as density functional theories p8|^3| and computer simulations |Tl|,^^-^ have 
the potential to investigate in detail the role of microscopic interactions and to resolve the 



fine structures of tlie fluids under consideration on a molecular scale. Whereas the former 
theories are based on severe simplifications and focus on idealized systems on a macroscopic 
scale (i.e., /im and larger), the latter require high computational efforts such that only rather 
small systems like pores and slits can be studied. Our present analysis is supposed to pro- 
vide a link between these two approaches in that it keeps track of the microscopic molecular 
interactions between the various constituents but leaves out fine structures on the molecular 
scale. This allows us to investigate liquidlike structures on a mesoscopic scale commensurate 
with the actual experimental system sizes. As basic ingredients for the molecular interac- 
tions we adopt Lennard- Jones, i.e., long-ranged interaction potentials. These potentials are 
not only applicable for rare gases but resemble also reliable effective potentials for molecules 
like small alkanes or alcohols [Q. 



In order to provide the information and terminology required for the subsequent con- 
siderations in Sec. |l| we recall some basic results for wetting of fiat and homogeneous sub- 
strates |^;S^ and extend them to the case of homogeneous substrates covered by a laterally 



homogeneous overlayer. In Sec. |IT| we present the theoretical basis for our analysis of 
wetting of chemically structured substrates from which we derive many analytical and nu- 
merical results in Sees. |^ and 0, respectively. These results extend and complete earlier 



studies which also focused on the mesoscopic structure of liquidlike layers on heterogeneous 
substrates [^,^]. Our analysis is carried out within mean field theory without taking into 
account capillary wavelike fluctuations and the strong fluctuations arising near critical points 
of the fluids. We focus our analysis on three paradigmatic chemical surface structures: a 
chemical step generated by two adjacent quarter spaces fllled with different materials, a 
chemical step generated by two different adjacent surface layers covering a homogeneous 
substrate, and a chemical stripe generated by a slab of different material immersed in an 
otherwise homogeneous substrate orthogonal to the common flat surface. Finally, in Sec. |V| 
we extend our approach to arbitrary chemically structured substrates and summarize our 
results in Sec. |Vli| . Technical details of the substrate potentials and of the density functional 
description are presented in Appendix ^ and p, respectively. 

II. WETTING OF LATERALLY HOMOGENEOUS SUBSTRATES 

A. Density functional theory 

Our theoretical analysis is based on the following simple |5^ but nonetheless success- 



ful [Q density functional theory for a spatially inhomogeneous number density p(r) of fluid 
particles: 

J^([p(r)];T,/x)= [d'rfHs{p{r),T)+ / rfV[V(r) - /i]p(r) 
Ja Ja 

+ - f d^r f dh-'w{\r-r'\)p{r)p{r'). (2.1) 

2 Ja J a 

The minimum of f2[p(r)] with respect to p(r) yields the grand canonical free energy of the 
fluid corresponding to a prescribed temperature T and chemical potential fi. V{r) denotes 
the substrate potential and w{r) is the attractive part of the fluid interparticle potential. 



A is the volume occupied by the fluid particles, fns is the free energy density of a hard- 
sphere reference fluid. This contribution takes into account the repulsive part of the fluid 
interparticle potential within a local density approximation which neglects short-ranged 
particle-particle correlations. These correlations become important in the close vicinity 
of the substrate surface. Nonetheless Eq. ( |2.1| ) has turned out to provide a very useful 
description of wetting phenomena in the case that the adsorbed liquidlike fllms are much 
thicker than the diameter o"/ of the fluid particles. 

The fluid particles are assumed to interact via a Lennard- Jones potential 



0/(0 = 4e/ 



r 



(2.2) 



We apply the Weeks- Chandler- Andersen (WCA) procedure [Q to split up this interaction 
into an attractive part (pattif^) and a repulsive part 0rep(^)- The latter gives rise to an 
effective, temperature dependent hard sphere diameter 

21/6,, 

(prepir) 



.(T)^ y *|l-exp^-rg^j| (2.3) 



which is inserted into the Carnahan-Starling approximation for the free energy density of 



the hard- sphere fluid ^A\ 



fHsip, T) = ksTp (inipX^) - 1 + ^Y3^) (2.4) 

where t] = ^p{d{T))^ is the dimensionless packing fraction and A is the thermal de Broglie 
wavelength. 

The attractive part of the interaction (f)att{f) is approximated by 

,I;(r) = ^^(r2 + a2)-3 (2.5) 



with 

Wo 



I dhw{r)= [ £r(l)auir) = -^V27iefa^r. (2.6) 

Although Eq. ( p.5|) does not strictly implement the WCA procedure corresponding to 
Eq. (|2.2|), it resembles it closely and offers the valuable advantage of reduced computa- 
tional efforts because the form oiw{r) allows one to carry out certain integrations over w{r) 
analytically. For large particle separations r one has w{r -^ oo) ~ r^^; the amplitude is 
chosen such that the integrated strength equals that of the attractive contribution (f)au ob- 
tained by the strict application of the WCA procedure. The third integral in Eq. ( p.l| ) takes 
into account the attractive fluid interparticle interaction within mean fleld theory based on 
Eq. 



B. Bulk phases 

In a bulk system the particle density p^ (where j = l,g denote the liquid and vapor 
phase, respectively) is spatially constant, leading to (see Eq. ( p.l|) ) 

^bip^, T, fi) = fHsip-y, T) + -wop^ - /up^ (2.7) 

for the grand canonical free energy density of a bulk fluid. Minimization of Vtb with respect 
to p^ yields the equilibrium densities. The conditions for liquid-vapor phase coexistence 
p = po(T) are 
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on,, 



P=P. ^P 



= and ^^{pg) = VL^{pi). (2.8) 

p=pi 
Off coexistence, i.e., for p ^ pq, only the liquid or the vapor phase is stable. In this case the 



dp 

, fo] 
density of the metastable phase corresponds to the second local minimum of fib. 



C. Wetting of laterally homogeneous substrates 

Here we recall some basic results for wetting phenomena (for reviews see, e.g. 



Refs. [50,55- pq] ) as they follow from Eq. (|2.1|) for a homogeneous substrate within the 
so-called sharp- kink approximation (see Fig. |1| and Ref. |Q). This allows us to formulate 
our corresponding findings for substrates which are heterogeneous either in the direction 
normal to the surface or in lateral directions. 

In the following we consider a flat substrate located in the half space ty = {r G M^l-z < 0}. 
The substrate is either homogeneous or composed of a homogeneous part wh = {r G M.^\z < 
—dgz} covered by a homogeneous surface layer ws = {r G M^| — dgz < z < 0} oi different 
chemical species. It consists of particles located on an orthorhombic lattice with lattice 
constants Qi in the i direction (i = x,y, z). For reasons of simplicity the lattice constants gt 
are assumed to be constant throughout the whole substrate. The surface layer consists of 
n = d + 1 monolayers with the uppermost and the lowest monolayer located at 2; = and 
z = —dgz, respectively. The theoretical description (Eq. ( p.l|) ) is identical for both substrate 
types; it only differs with respect to the considered substrate potential V{z). The substrate 
acts as an inert spectator phase. The interaction potential between a fluid particle and an 
individual substrate particle is taken to be of the Lennard- Jones type, too: 



^H,s{r) = 4:eH,s 



ch,s\^'^ (<^h,s 



(2.9) 



where H and S denote the molecules in wh and ws, respectively. The substrate potential 
V{r) = V{z) = Vau{z) + Vrep{z) follows as a laterally averaged pairwise sum of (J)h,s over all 
substrate particles, resulting in 
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wf uf 



att\ 



{z + id+l)gzr iz + id+l)gzr 



{z + dgzYJ " \z^ {z + dg 
2{d + l)ul - {2d + l)u^ 



< 

z^ z^ 



+ 0(;z-^), z:^dgz (2.10) 



for the attractive potential. In Eq. ( p.lO|) we have used the relations M4 ' = f^/z^s ' ■ 
Moreover, the coefficients can be expressed in terms of the parameters eH,s and (yH,s of the 
molecular interactions and the lattice spacing g^. The repulsive contribution is 
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+ 0(z-i"), z:^dg,. (2.11) 



In the case of a homogeneous substrate, uf = uf = Ui, V{z) reduces to 

n^) = -E|- (2-12) 

i>3 

The limit (i — i> 00 corresponds to a homogeneous substrate composed of the chemical species 
of the surface layer. 

For the above model the full minimization of the functional in Eq. (|2.1| ) with respect 
to the spatially inhomogeneous, smooth density p{z) can be carried out only numerically. 
However, it has turned out that the minimization restricted to the subspace of piecewise con- 
stant density profiles provides a surprisingly accurate analytic account of the corresponding 
effective interface potential ||59|j60|]. Within this so-called sharp-kink approximation (Fig. fj). 



p{z) = Q{z - d^)[Q{l - z)pi + Q{z - l)p,l (2.13) 

pi and Pg are the number densities of the bulk liquid and vapor phase, respectively (see 
Subsec. [II B|) . For p < po pg is taken to be the actual vapor density off coexistence and pi 
is the density of the metastable liquid phase. 9 denotes the Heaviside step function, and 
/ is the thickness of the adsorbed liquidlike wetting layer, i.e., ^ = / is the position of the 
emerging liquid-vapor interface. In terms of the actual smooth density distribution p{z) the 
position of the interface can be defined by, e.g., p{z = I) = ^{pi + pg). The length d^, takes 
into account the excluded volume for the centers of the fluid particles due to the repulsive 
part of the substrate potential. We assume that d^ = ^{a^ + CTf). The insertion of Eq. ( p.l3| ) 
into Eq. ( p.l|) yields 



nmz)]-,T,p-,[w],[V])=Anb{pb-,T,p) + An,{[p{z)]-T,p-[w],[V]), (2.14) 

i.e., a decomposition of the grand canonical functional into bulk and surface contributions. 
A is the volume occupied by fluid particles. In the thermodynamic limit this is the half space 
{r = (ry, z) G R^\z > 0}. A is the area of the substrate surface, flf, is given by Eq. ( p.7|) 
whereas Qg yields the so-called effective interface potential fis(/): 

nsMz)]) = ^s{l) = AQtl + cr^i + crig + uj{l). (2.15) 



The flrst term in Eq. ( p.l5| ) measures the cost in free energy for forming a liquidlike wetting 



layer of thickness / if in the bulk the vapor is the stable phase. One has AVt^ = ApAp + 
0{{Ap,)^) where Ap = po - p. With 

6 



t{z) = dz' d^r\\w{Jrl + z'^) (2.16) 

2 

as the potential of a fluid particle interacting with a half space filled with the same fluid 
particles, within the above sharp-kink approximation one has for the substrate-liquid surface 
tension 



oo 





and 



<^wi = -\p1 j dzt{z)+pi j dzV{z)-d^n^^ (2.17) 



oo 

aig = ~{Apf fdztiz) (2.18) 



for the liquid-vapor surface tension. The /-dependent part uj{l) of the effective interface 
potential is given by 



uj{l) = Ap{pi / dzt{z)- I dzV{z)) . (2.19) 




Minimization of the interfacial free energy VLs(l) with respect to / yields the equilibrium 
film thickness /q and the substrate- vapor surface tension a^g = minfis(0 = ^s(/o)- From 

Young's equation one obtains for the contact angle cos(6') = {a^g — (Jwi)l<yig = 1 + ^{!'o)/o'ig 

which is thermodynamically well-defined only at two-phase coexistence Ap = 0. 

From Eqs. ( p.5|) and (|2.16| ) one has t{z) = — ^tj z~^ and therefore 

i>3 

^(0 = E|- (2-20) 

(Here we do not consider terms ~ Z~^ln/ generated by van der Waals tails in the density 
profiles [60(c)] not captured by the ansatz in Eq. ( p.l3| ).) The coefficients of the leading 
terms are 02 = -^{u^ — tsPi), which is known as the Hamaker constant, and 03 = -^(2((i -(- 
l)uf — {2d + l)uf — (t4 + ?>tzdw)pi)] for a homogeneous substrate with Uj=uf = Uj one 
has 03 = -^(^4 — (^4 + 3t3du))pi). We note that the Hamaker constant for a homogeneous 
substrate covered by a homogeneous surface layer does not depend on the properties of the 
surface layer but all subdominant terms do. 

The substrate is said to be completely wet by the liquid phase if, at coexistence Ap = 0, 
Iq is infinite. Since a;(/o = cxd) = one has a^g = cTwI + o'lg and 6' = for completely 
and (Tuig < (Twi + (Tig and 6^ > for partially wet substrates. If at T = T^ the thickness 
lo(T,p = po) jumps from a finite value for T < T^; to /q = 00 for T > T^ the system 
undergoes a first order wetting transition. If the film thickness grows continuously upon 
approaching the wetting temperature, i.e., 1{T — > Tw,p = po) —^ cxd the system exhibits a 



critical wetting transition. The necessary condition for critical wetting is that 02 (T) changes 
sign at Tw from a2(T < T^) < to a2(T > Tyj) > and a3(T = T^) > 0, i.e., M4 — 
(^4 + ?>t^dyj) piiTyj) > 0. The critical wetting transition temperature is given implicitly by 
u^ = t3pi{Tw), or equivalently M3 = t^piiT^) for a homogeneous substrate. Even with 
an additional surface layer, which modifies the substrate potential in the vicinity of the 
substrate surface, the transition temperature of a critical wetting transition is determined 
only by the properties of the underlying homogeneous substrate wh- This implies that two 
substrates which differ only with respect to their overlayers have the same wetting transition 
temperature if the wetting transition is continuous. From the above formulae one can infer 
that the occurrence of a critical wetting transition on a planar and homogeneous substrate 
hinges on the subdominant contribution ~ z''^ in the asymptotic expansion of V{z) for 
large z because for the attractive fluid-fluid interaction as given by Eq. ( |2.5| ) ^4 = and 
ta = — ^ifoCj > 0. It is possible to fulfil the necessary condition a^{Tw) > for critical 
wetting by choosing an appropriate surface layer. In the case of critical wetting the film 
thickness diverges as lo{T y T^) = —803/202 ~ (Tw — T)~^. Irrespective of the order of the 
wetting transition, upon approaching coexistence along a complete wetting isotherm one has 
lo{T > T^, A/i \ 0) = (2a2/A(]6)i/3 _ (A/i)-i/3. For T > T^ the Hamaker constant 02 is 
always positive. The leading divergence of Iq for a complete wetting transition is independent 
of any different surface layer covering the substrate. 

III. MODELS FOR WETTING OF STRUCTURED SUBSTRATES 

A. Simple chemical step 

As a basic element for more complicated structures we first analyze the wetting properties 
of a substrate which exhibits a single, simple chemical step (SCS), i.e., a flat substrate 
composed of two adjacent quarter spaces filled with different chemical species (see Fig. |^ 
with n = cxd). The substrate particles occupy the half space w = {r G M'^|z < 0} and the 
heterogeneity defines the position x = 0. The system is translationally invariant along the 
y direction. The substrate and the fluid particles interact via Lennard- Jones potentials 



^±{r) =4ed 



r / 



(3.1) 



where the "+" and "— " signs refer to substrate particles located in the quarter spaces 
U74. = {r G R'^|x > A 2 < 0} and W- = {r G R^|x < A 2 < 0}, respectively. The substrate 
potential V{x,z) = Vatt(yX,z) + Vrep{x,z), as obtained by a pairwise summation over all 
substrate particles, is given by Eqs. ( |A^) and ( [A^ ) in Appendix 0. In the limit |a:| — > cxo, 
V{x, z) asymptotically approaches the substrate potential of the respective homogeneous 
substrate: 

l/(. ^ ±«3. ^) = -f - I + I + ^^ + 0(X-') 

= '4(.) + ^^^ + 0(x-') (3.2) 



where the coefficients of the subdominant terms omitted here depend on z. Thus the ad- 
sorbed hquidhke wetting film exhibits the asymptotic thicknesses /+ for x ^ oo and /_ 
for X -^ — oo which are determined by the properties of the homogeneous substrates "+" 
and "— ", respectively. In general one has to take into consideration two different wetting 
temperatures T^ and T~ , one for each of the two semi- infinite, homogeneous substrates. 

With the substrate potential F(a;, z) given, in principle the functional in Eq. ( p.l|) can 
be minimized with respect to the full density distribution p{x,z). However, this can only 
be done numerically and requires a huge computational effort. Therefore we focus on low 
temperatures which allows us to restrict the minimization to steplike proffies 

p{x, z) = [Q{-x)Q{z - rf^) + e{x)e{z - d+)] 

X [p,e(/(x) -z)+ pgQ{z - l{x))]. (3.3) 

The quantities d^ = \{<J± + erf) take into account the different excluded volumes in the 
vicinity of the substrate surface. Because of the rapid decay of the repulsive forces d^ is 
taken to vary steplike at a; = between c?+ and d~. The approximation used in Eq. ( p.3|) 
does not capture the fluid density oscillations very close to the substrate surface. However, 
these fine structures of the particle density are expected to have only small effects if the 
liquidlike wetting ffims are rather thick, as it is the case in the vicinity of a critical or a 
complete wetting transition. The sharp-kink approximation yields an exact prediction for 
the transition temperature Ty^ of a critical wetting transition which remains unchanged if 
more sophisticated models are applied [^. Although within the sharp-kink approximation 
the predicted thickness of adsorbed liquidlike wetting layers is not quantitatively accurate 
this approach is expected to capture the essential features of the wetting phenomena as 
considered on the present mesoscopic scale. 

The grand canonical free energy functional - which via Eq. (|3.3| ) is a functional of the 
function l{x) that describes the local position of the liquid- vapor interface - can be system- 
atically decomposed into bulk, surface and line contributions: 

ni[p{x,z)];T,p;[w],[V])=An,{pg,T,p) + An,il^;T,p;[w],[V]) 

+L^Qi{[l{x)];T,p-[w],[V]). (3.4) 



The explicit expressions for these contributions are given in Appendix ^. In Eq. 
A = LxLyLz is the volume filled with fluid particles, A = L^Ly is the surface area of the 
substrate surface, and Ly is the linear extension of the chemical step. Qb, given by Eq. (p. 7]), 
is the bulk free energy density corresponding to the stable bulk vapor phase. The surface 
contribution 

fi,(/±) = i(l]+(/+) + fi;(/_)) (3.5) 

is the arithmetic mean of the surface free energy densities corresponding to the substrates w+ 
and W- covered by liquidlike films of thickness /+ and /_ which are exposed to the potentials 
V+{z) and V-{z), respectively (see Eq. (|2.15|) ). l± minimizes the corresponding Eq. (|2.15|) 
with pg and pi as the bulk vapor and the metastable bulk liquid phase, respectively. 
The line contribution Qi is due to the substrate heterogeneity and reads 

ni[i{x)] = T{dti±) + m^)]- (3.6) 



The expression T{d^,l±) does not depend on the profile l{x). Minimization of the nonlocal 
functional (D[/(x)] (which is given in Eq. ( p^) ) yields the equilibrium liquid- vapor interface 
profile l{x) and the line tension 

T{T,fi) = min ni[l{x)] = ni[I{x)]. (3.7) 

{l(x)} 

A gradient expansion of (I;[/(x)] in leading order leads to the local functional uJioc[l{x)] (see 
Eq. (|B9|) ) that provides a prescription of how to express the functional expressions given 
by a simple, phenomenological interface displacement model in terms of the microscopic 
parameters of the underlying molecular interactions. 

The Euler-Lagrange equation (ELE) following from the functional derivative of the non- 
local functional a;[/(a;)] is 



5ni[i{x)] 



5l{x) 



_ 5Co[l{x)] 

1=1 ^Kx) 



1=1 



Afifc - Ap[p, t{l{x) - rf+) - V{x, l{x))] + /(x, l{x)) 
+Appi {i{x, T{x) - d+) - i{x, l{x) -d-)} = (3.8) 



with t{x, z) given by Eq. (|B7|) and 



oo l{x')-l{x) 

I{x,l{x)) = {Apy dx' I dz' w{x - x' , z') (3.9) 

-oo 



where w{x, z) is given by Eq. ([B8|) . The ELE is a nonlocal integral equation for the function 
l{x). Within the local theory the double- integral in Eq. ( p.9| ) is replaced by a differential 
expression leading to 

+Appi {t(x, l{x) - dl) - i{x, l{x) - C)} ■ (3.10) 

Equation ( |3.10| ) is often referred to as "augmented Young equation" ||HI[ . 



B. Chemical step within a surface layer 

The analysis in the previous subsection requires that the substrate is composed of two 
adjacent halves which themselves are homogeneous along both the x and the y direction. 
However, it is not necessary that the substrate halves are chemically homogeneous in the z 
direction. 

This allows us to consider within the same formalism a homogeneous substrate wh 
covered by a surface layer which itself is composed of two adjacent homogeneous layers 
ws,± = {r G M^|a; ^ A —dgz < z < 0} of different chemical species. In the following 
we refer to this type of substrate as a feyer chemical step (LCS). The two halves of the 
layer meet at x = (see Fig. |^). Both surface layers consist of ra = rf + 1 monolayers 
with the uppermost and the lowest monolayer located at z = and z = —dg^, respectively. 
This model mimics rather closely the kind of substrate inhomogeneities generated by, e.g., 

10 



microcontact printing (see Sec. 1). The expression for the substrate potential V{x,z) of 
this system is rather comphcated (see Eqs. (|A4|) and (^)). As mentioned in Sec. P^IC 



the leading contribution to the attractive substrate potential is only determined by the 
underlying homogeneous substrate wh- Thus in the case that both substrate halves undergo 
a critical wetting transition there is only a single common wetting temperature T^ for the 
whole substrate because the Hamaker constant is only determined by wh- 

Far from the heterogeneity, for fixed z, the potential V{x, z) asymptotically approaches 
the substrate potentials V±{z) of the two semi- infinite, laterally homogeneous substrate 
halves covered by a surface layer via 



3{d + l) {u2 -u^) 1 ^ ^^^^-5^ 
8 x^ 



V{x -^ ±oo, z) = V±{z) + sign(x) ^ ^;^ ^^^ + C(x-^). (3.11) 



Again the coefficients of the subdominant terms omitted here depend on z. From Eq. ( p.ll| ) 



it follows that the wetting layer thickness l{x) asymptotically approaches the constant values 
l± corresponding to the respective decorated substrate. 

C. Chemical stripe 

As a third system we consider a substrate exhibiting a chemical sfripe (CST). The stripe 
is realized by insertion of a slab Wgt = {r G M^||x| < a/2 A 2; < 0} of different chemical 
species "+" into an otherwise homogeneous substrate w = {r e M^l^ < 0} composed of 
particles denoted as "— " such that in top view a chemical stripe is formed (see Fig. ^. 
The system is again translationally invariant in the y direction. The substrate potential 
V{—x, z) = V{x, z) is given by Eqs. ( [A^ ) and ( [AT]) . For z fixed in the limit of large \x\ one 
has 

V{\x\^oo,z) = -%-^ + %-^ <eZ^ + o(,-5) (3.12) 

Z-^ Z^ Z^ Qx X* 

(where again only the coefficients of the subdominant terms depend on z) implying that the 
equilibrium liquid-vapor interface is also symmetric with respect to x = 0, i.e., l{—x) = l{x), 
and that it asymptotically approaches the constant value Z_ for |a;| ^ 00 determined by the 
properties of the homogeneous, flat substrate w. For a wide stripe width one has 



11 ?/ ?/ 

V{x,z) = |-^ + -^ + 0(a-='), a > (7, and |a;| < a/2, (3.13) 

Z"^ z^ z^ 



so that for small |x| the proflle l{x) is close to the fllm thickness /+ of a liquidlike fllm on 
the homogeneous and planar "+" substrate. 

The sharp-kink density profile used for the analysis of the CST is 

p{x, z) = {0 (|x| - ^) e{z - rf-) + (I - |x|) e{z - di)} 

X [0(/(x) - z)pi + e{z - l{x))pg] . (3.14) 

Inserting p{x, z) in Eq. ( p.l|) yields 
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fi([p(x,z)];T,/i; [w], [V]) = An,ipg,T, fx) + Afi,(L;T,/i; [w], [V]) 

+Lyni{[l{x)];T,fi;[w],[V]). (3.15) 

The surface contribution Qg (see Appendix P) is determined by the properties of the ho- 
mogeneous, flat substrate w. Thus in contrast to the SCS here one only has to deal with 
one wetting transition on w at T = T~ . The wetting transition on the stripe part Wgt is 
suppressed due to the finite lateral extension of the stripe. 

The effect of the stripe on the liquid- vapor interface is captured by the line contribution 
to the free energy functional. This contribution reads 

ni[iix)] = T{dti~) + mx)]- (3.16) 

Whereas T{d^,l_) (see Eq. (|B14| )) is independent of /(x), tD[/(x)] (see Eq. ( pi5|) ) is a func- 
tional of l{x). Functional differentiation of i7;[/(x)], i.e., of cj[/(x)], yields the ELE 

a/2 l(x)—d^ 

AH, - Apip, tm - 4) - V(.. '»)1 - AP P, / i^' j i^ ^(^ - < A 

-a/2 l{x)-d+ 

oo T{x')-T{x) 

+ {Apf I dx' I dz'w{x-x',z') = 0. (3.17) 

-oo 

The corresponding ELE within the local approximation uJioc[l{x)] of cj[/(x)] is 

An, - Ap[pi t{l{x) - 4) - V{x, l{x))] 

a/2 T{x)—d^ _ 

-App/j dx' j ^--(---^-)-(^f|!j^ = o• (3.18) 

-a/2 l{x)-d+ 

In the following we always discuss the actual equilibrium liquid-vapor interface profiles; 
therefore we omit the overbar. 



IV. PROPERTIES OF THE INTERFACIAL PROFILES 

A. Curvature behavior 

Within the local theory the ELEs for the SCS, the LCS, and the CST determine the 
local curvature K{x) of the planar curve {x, l{x)): 

"'^"^("^ = (1 +7/(1^)^)3/^ = ^^"' ^^")) ^'-'^ 



with 



f/(x, z) = An, + ^^(''l^'J^t^ + c{x, z- di). (4.2) 
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The "local" effective interface potential u!{x,z;d^) is given by Eq. (p6[). The function 
C{x, z; d^) takes into account the difference between d^ and d~, i.e., it vanishes for (i+ = d~ 
and for (i+ 7^ d~ in the limit |x| — >• cx) (see Eqs. ( p.lO| ) and ( p.l8| )). In the limit |a:| — > cx) the 



lines z{x) defined implicitly by U{x, z) = asymptotically approach those values for which 
the function AQbl + uj±{l) (compare Eq. (|B2[)) is extremal with respect to /. One of these 
values corresponds to the global minimum and thus is the equilibrium film thickness l±. In 
addition to this one there may be more extremal values depending on whether the system 
is at or off coexistence and on the type of wetting transition under consideration. Since l{x) 
asymptotically approaches the values l± for |a;| — > cxd it also approaches the contour line 
given by U{x, z) = corresponding to the global minimum of AQf, I + uj±{l). Figure ^ shows 
an example for the line U{x,z) = on a SCS and the corresponding interfacial profile l{x) 
for a thermodynamic state along a complete wetting isotherm. The sign of the curvature is 
given by the sign of U{x, z). The point where l{x) and the line U{x, z) = intersect is the 
turning point of the profile. Due to this curvature behavior the function l{x) for a SCS or 
a LCS is monotonous. Analogously, for a CST l{x) is monotonous in each of the intervals 
X < and a; > 0. 

In general the curvature behavior as predicted by the actual nonlocal theory is differ- 



ent from that obtained within the local approximation. However, in Ref. |B^ we have 
demonstrated that the differences between the local and the nonlocal results are negligible 
if the local curvature is small, as it is the case for complete and critical wetting transitions. 
Therefore, although the behavior of the integral expression in Eqs. ( |3.8| ) and (|3.17|) is not 
analytically transparent, we expect the curvature behavior of the profiles upon approaching 
a critical or complete wetting transition to be practically the same within the nonlocal and 
the local theory. 

B. Asymptotic behavior of the interface profiles 

The asymptotic form of the substrate potential V{\x\ -^ 00, z) of a heterogeneous sub- 
strate gives rise to a behavior 

/(|x| ^ 00) = /± + 6l{x) with 6l{x) = X x-" (4.3) 

of the interface profile with a characteristic exponent a. A is the amplitude of the so-called 
van der Waals tails (VDWT) 6l{x) which can be determined analytically for all three sub- 
strate types SCS, LCS, and CST for those cases in which parts of the substrate, considered 
per se in the corresponding homogeneous limit, undergo a critical (T y T^,Afi = 0) or a 
complete (T > T^, Afi \ 0) wetting transition. These results are based on the assumption 
that the system is sufficiently close to the wetting transition temperature or to two-phase 
coexistence such that l{x) ^ d^ and thus 



Al],+ ^"^^'^) 



= I^ioc){x,l{x)) (4.4) 

l=l{x) 



dl 

for all three substrate types. uj{x,l) is given by Eq. ( |B^ ) omitting the arguments d^. 

Since with l{x) ^ c?^ also l± ^ d^ Eq. ( [4.4| ) can be used for inserting the ansatz 
l{x) = l± + 6l{x) into its local version. Considering then the behavior of Eq. i\i.4\ ) for large 
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|x| and exploiting the fact that l±, which corresponds to the hmit |x| -^ oo, minimizes 
Eq. ( [BID leads to the following expansion: 

a5l{x) + bx~'^ = aig5r'{x), a > 0, (4.5) 

with coefficients a and b. The term ~ x~" stems from the leading term in V{\x\ -^ oo)—V±{z) 
(see Eqs. ( p.2|) , (p.ll|) , and ( p.l2| )) which determines the value of a. The leading asymptotic 



behavior of the solution of Eq. ( [4.5| ) is given by 

Sl(x) = --x-''. (4.6) 

a 

Within the nonlocal theory the right-hand side is replaced by the leading order in the series 
expansion of the integral I{x,l± + 61 (x)): 

a6l{x) +bx~'' = I{x,l± + 6l{x)) ^ x-^, (3 > 0. (4.7) 

Thus the asymptotic solution 6l{x) in Eq. (|4.7|) for the VDWT is equal to that obtained 
within the local theory (Eq. ( |4.6D ) if the right-hand side in Eq. ( |4.7| ) is subdominant as 
compared with the left-hand side, i.e., ii (3 > a. For P = a the amplitude A (see Eq. (^^)) 
obtained from the local theory differs from that following from the nonlocal one. It turns 
out that this is the case for the LCS. If one would have f3 < a which, however, is not the 
case here, even the exponent of the power-law decay of 6l{x) as obtained within the nonlocal 
theory would be different from that within the local theory, i.e., the VDWT would also differ 
qualitatively. 

With the definition l{x -^ ±oo) = /± =F Sl±{x) for the SCS it follows that 

6h{x) = Xt,,x-' with Xfru = ±i^^^ (4.8) 

in the case of a critical wetting transition on the semi-infinite, homogeneous substrate w^ 
or w_. Along a complete wetting isotherm the result is the same up to a numerical factor: 

6l±{x) = Xi^pX-^ with Xi^p = ±f^^^^. (4.9) 

Both Eq. ( [4.8| ) and Eq. ( |4.9| ) are valid within the local as well as the nonlocal theory in 
agreement with Ref. |^ . For the LCS one finds a different power law which reflects the fact 
that for this type of substrate the lateral dependence of the potential enters only into the 
subdominant terms: 



Sl{x) = Xx-^ (4.10) 



with 



(a+-a3-)(l-|t3^) (nonlocal) 



Ail 

32 \a2\ 



At.t = r""' . /i ^; ,^ ^"'"'^ \. ., ' (4-11) 



^ * 'al — a^ ) (local) 



and 
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Xi^^ = { ^-^""^^ "/^^ (^ - 3^3^) (nonlocal) _ ^^_^^^ 

i^(4-«3) (local) 

Both for a complete and a critical wetting transition the local and the nonlocal theory yield 
the same power-law behavior but different amplitudes. Analogously, for the CST one finds 

l{\x\ — >• oo) = /_ + Sl{x) where 

3a ^4 4-%" (4_^3) 



Sl{x) = XcritX ^ 


with 


Xcrit — 


critical wetting and 






0l[X) = AcompX 


with 


^comp 



2g: 



,4 Cb'i flq 



/4 ^ _ ^ (4.14) 

2g^ a^ 

for a complete wetting transition on the substrate w. Thus for the LCS the VDWT do 
resolve the difference between the local and the nonlocal theory, but not for the SCS and 
the CST. 

On a laterally homogeneous and flat substrate the film thickness is given by l± = 
Sog /2|a^| for critical wetting and by l± ~ (a^/A/i)^/^ for complete wetting. For the deriva- 
tion of the expressions for the VDWT we have used these relations which do not hold for 
a^ = in the case of critical wetting and a^ = in the case of complete wetting, respec- 
tively. In these cases of complete wetting at Ty^ or of tricritical wetting higher-order terms 
have to be taken into account in order to determine the amplitudes of the VDWT. Finally 
we note that due to the behavior of the curvature of l{x) - as stated in Subsec. IV A - the 
signs of the amplitudes of the VDWT are fixed by sign(/4. — /_). 

C. Partial versus complete wetting of a chemical step 

The SCS and LCS allow one to realize thermodynamic states for which one half of the 
substrate, say w+, is completely wet whereas the other half is only partially wet. Such a 
state is realized for T+ < Tq < T~ and A/i = 0. The corresponding interface profile does 
not depend on whether this state is reached along coexistence (T /" Tq, A/i = 0) or along a 
complete wetting isotherm (T = Tq, A/i \ 0) (see Fig. ^. The liquid- vapor interface attains 
a finite value I- for x ^ — 00 and diverges for x — > 00. 

The characteristic length scale for lateral variations of the interfacial profile l{x) is the 
lateral height-height correlation length C,\\- For complete wetting of the substrate w± this 
length is given by (compare Ref. [ |63|] ) 



^ II. comp V '9 



^2 ^ -1/2 



,comp V 'y \ Jj 



l=l±/ V ^""2 

whereas for a critical wetting transition 



-^/|~(A/.)-2/3 (4.15) 

bar, 



tcrU=J^lT-t-J'^ (4.16) 



Sag 
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where t^ = {Tw—T)/Tw (compare Ref. 0]). In the case that the substrate half w+ undergoes 



complete and the substrate half w_ incomplete wetting we make the scaling ansatz 

/(x, A/x) = /+(A/i)/(x/e+,_^(A/i)) (4.17) 

for the evolution of the profile l{x, A/i > 0) from a shape which attains the constant value 
/+ in the limit x — > cxd for A/i > to the shape for A/x = with l{x — >• cxo, A/i = 0) = cxo. 
With this ansatz and in the case that the slope of the profile is small such that the left-hand 
side of Eq. ( |3.1CI| ) can be replaced by its lowest-order expansion aigl"{x) one obtains the 
differential equation 

i-f{yr^^~'^ = iC-^)ny) (4.18) 

for the scaling function f{y) in the case that the attractive fluid- fluid and substrate-fiuid 
interactions decay as r~^; f{y -^ oo) = 1 independent oi(. This demonstrates that the above 
scaling ansatz does indeed hold. In the limit A/i \ which corresponds to ^i^'^g^ — *■ oo, 
i.e., 1/ — ^ 0, from the asymptotic behavior of the solution of Eq. ( |4.18| ) it follows that 



/(a;^oo)=7co^^,a;2/('^+i) with ^.^^ = ("i^lilM") . (4.19) 

The coefficients in Eq. (|4.19| ) are defined by the asymptotic behavior of the effective interface 
potential iu+{l ^ oo) = a'^ l~°' + a+ /"" + C(/"''"^) where a = ( - A and k = ( - 3. For 
Lennard- Jones potentials cr = 2 and k = 3 so that 

/(x ^ oo) = (^^^ x^/^ A/i = 0, T+<T<TJ. (4.20) 



This power law has been predicted originally by de Gennes |S^, the expression for the 
amplitude is in accordance with Ref. [0. Figure |] shows the scaling function f{y) for 
C = 6. 

For the case T y Tq = T^ with the critical wetting transition temperature T+ on the 
substrate half w^ and AyU = the scaling ansatz 

l{x ^ oo,T) = /+(T)(7(x/e+,„,(T)) (4.21) 

leads to the differential equation 

9{y)-^''~'^-9{yr^''-'^ = g"{y) (4.22) 

for the scaling function g[y) with giy —>■ oo) = 1 independent of (. In the limit T y T+ 
,^|t ^^(T) diverges and from the behavior of the solution of Eq. ( [4.22| ) for y ^ one finds 

/(x^+oo) = 7c„ix2/(''+2) with 7cru= r \J^ M (4.23) 

for T = T+ so that for ^ = 6 
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l(x^oo)=^-j^j x^l\ A/i = 0,T = T+<r-. (4.24) 

Figure ^ shows the scahng function g[y) for C = 6. 

All these equations hold both for the SCS with a^ = -i^Wt ~ ^3Pi) ^^^ a'^ = -j-i'^t ~ 
{U + ?>hd^)pi) and for the LCS with a+ = ^(wf - ^pi) and a^ = ^{2{d + 1)m| - {2d + 
l)u^ — {ti+'it^dyj) pi) . The scaling functions / and g are independent of the thermodynamical 
parameters T and p and of the amplitudes of the molecular interactions. These dependences 
are completely absorbed into ^t and /+. However, / and g do depend on the exponent C 
of the power-law decay of the attractive interparticle potentials. We note that the scaling 
functions / and g do not take into account the lateral long-ranged VDWT of l{x) (see 
Subsec. [IV B|) . These VDWT give rise to additional small corrections to the behavior of 
fiy -^ oo) and g{y -> oo). 

D. Partial versus complete wetting of a chemical stripe 

A similar analytic calculation as in the previous subsection can be carried out for the 
case of a CST with the parameters chosen such that at coexistence Ap = and a fixed 
temperature T+ < T < T~ the outer area is only partially wet whereas the stripe region is 
completely wet in the limit a = oo. It turns out that the analytical solution for the shape 
of the profile in the case of large a is a semi-ellipse: 

i.e., the half axes are r^ = a/2 and r^ = (2a^/o"/(,)^/^A/a/2. The corresponding excess 
coverage, in comparison with the case a = 0, scales as T^xicb -^ C)o) ~ a^^'^. In the limit 
a ^ oo at both stripe boundaries x = ±a/2 one recovers the square-root behavior given by 
Eq. ( [4.20| ). In this limit the presence of the second, distant boundary of the stripe gives rise 
to corrections to this square-root behavior of the interfacial profile: 



l{x; a -^ cxd) = I ) y/x 

(yig 



1/4 

(4.26) 



2 



1 - — - ^^^ + 0{{x/a 



2a 8a2 



The expression in Eq. ( [4.26| ) corresponds to a shifted coordinate system in which the two 
boundaries of the stripe are located at x = and x = a — > oo. 

E. Retardation 

Due to retardation for large separations r dispersion interaction potentials decay as r~^ 



rather than r~ |^8[. These retardation effects become important for thick wetting layers. 
Therefore they have to be taken into account in the discussion of the VDWT and of the film 
morphology for the crossover between partially and completely wetted substrate parts. For 
homogeneous and planar substrates in the retarded regime the substrate potential decays as 
V{z ^ Of) = —Vi z~'^ — f5 z~^ + 0{z~^\ If the interaction potential of the fluid particles is 
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also retarded one has t{z ^ cr/) = —S4Z "^ — s^z ^ + 0{z ^) so that the effective interface 
potential turns into 

uj{l » a;) = ^ + ^ + Oin (4.27) 

with 63 = -g^(t'4 — S4P1) and 64 = ^(^5 — (ss + 4s4C?^)p;). Thus for critical wetting 1{T y 
TJ = -454/353 with 63 changing sign at T = T^ and /(A/i \ 0) = (3 63/Afife)i/4 ~ 
(A/i)"^/*^ along a complete wetting isotherm. 

The substrate potential of a SCS, V"(a;, 2;) = Vatt{x,z) + Vrep{x,z) with Vatt{x,z) in 
accordance with Eq. (|A8|) , approaches 



z^ ■ 2x^ 



Vau{x - ±00, ;^) = -^ + -^.^ + 0(x-^) (4.28) 



for large \x\. Equation ( |4.28|) implies 5l±{x) = l{x —>■ ±00) — /-^ ~ x~'^ for the VDWT if the 
substrate w± undergoes a critical or complete wetting transition. Analogously, for the LCS 
and the CST one finds 6l±{x) = l{x — >• ±cxd) — /_|_ ~ x~^. 
For C = 7 Eqs. (|1|) and (^^ imply 

l{x ^ cx)) = f^l a;2/^ (4.29) 

in the case of complete wetting of w+ and partial wetting of W- (i.e., T^ < Tq < T~ , 
AyU = 0) and 



18 6^ 



1/6 
l{x ^ 00) = ( ^^^^ ) x'/^ (4.30) 



at critical wetting (i.e., T+ = Tq < T^ , A/i = 0). 



V. NUMERICAL ANALYSIS OF THE INTERFACE MORPHOLOGY 

Although several properties of the liquid- vapor interfaces can be determined analytically, 
even within the sharp-kink approximation the behavior of l{x) over the whole range of 
X values and for arbitrary values of T and fi can only be obtained numerically. Within 
the local theory for the SCS and the LCS we solve the ELE given by Eq. ( |3.1CI| ) and for 
the CST Eq. ( |3.18|) . The corresponding analysis within the nonlocal theory requires to 
minimize numerically the functional expression a}[/(x)] (see Eq. ( |B4D ) instead of solving the 
corresponding ELE. The reason is that as pointed out in Ref. [^ the numerical procedure for 



solving this nonlocal ELE leads to severe difficulties. In Ref. [^] the numerical techniques for 
calculating liquid- vapor interface profiles and a comparison of the results as obtained within 
the nonlocal and the local theory are described in detail. It turns out that in general these 
differences are very small. Therefore in the following we focus on the numerical analysis of 
the local theory, keeping in mind that these results are practically indistinguishable from 
those following from the nonlocal theory although, in principle, the latter is the more reliable 
one. 



A. Chemical step 

Figure ^ displays profiles of liquid-vapor interfaces for a substrate with a chemical step 
undergoing complete wetting of both w+ and W-. The parameters of the interactions are 
chosen such that both substrates exhibit critical wetting transitions, Wj^ at T^ = ksT^jef = 
1.0 and w_ at T* = 1.2. The profiles of the liquid-vapor interfaces are calculated for 
T* = 1.3 > T^* and different values of A/i. In the limit A/x -^ both asymptotic film 
thicknesses l± diverge according to l± ~ (A/i)~^'^. The interface profile is extremely broad; 
we note that the scales of the two axes differ by about three orders of magnitude. 

The crossover of l{x) from /_ to /+ is governed by two characteristic length scales, one 
for each substrate half w^ and w_. Whereas along the isotherm the position |a;o| ^ erf of the 
turning point remains practically unchanged, the two points x^ and x_ where the profile 
l{x) deviates from its asymptotes by 10% of |/+ — /_| diverge according to a power law. The 
characteristic length scale of the lateral variation of the profile is governed by the lateral 
interfacial height-height correlation length C,\\ given by Eq. ( |4.15| ). Figure ^ demonstrates 
that x± and ^jfcomp ^^^ proportional to each other. 

Figure P shows liquid- vapor interface profiles calculated for a SCS which exhibits a critical 
wetting transition at T^ = 1.0 for both substrates w_ and U7_|_. The interface profiles have 
been obtained for different T along the thermodynamic path {T —>■ T^j, Afi = 0). As before 
the position of the turning points remains fixed at |a;o| < o"/, the values of x+ and X- 
diverge for T — > T^, i.e., l± -^ oo. Here the divergence of x± is determined by Cf^^^^ given 
by Eq. ( |4.16| ). From Fig. |10|it can be seen that also in the case of a critical wetting transition 



x± and ,^11 ^^jj are proportional to each other. 

The behavior of /(|x| — > oo) is ultimately determined by the VDWT given in Eqs. 



and ( ^4.9| ). However, it turns out that these power-law tails in l{x) become relevant only for 
such large values of |x| where in fact 6l{x) is smaller than the numerical error. Therefore 
the analytically known VDWT are numerically not accessible. 

As demonstrated in Ref. [^ the more general nonlocal theory can be approximated by 
the corresponding local theory without losing significant quantitative accuracy. Since the 
slope dl/dx is small one can go even a step further and replace within the local theory the 
expression a/1 + {dl/dxY — 1 by {dl/dxY /2 to be used in the local functional uoioc[l{x)\ (see 
Eqs. ( |B4| ) and ( |B9| ) ) . This gives rise to the ELE of a square-gradient theory: 

aigl"{x) = AQ, - Ap[pit{l{x) - d+) - V{x,lix))] 

+Appi {i{x, l{x) - d+) - i{x, l{x) - d~)} (5.1) 

which is the analogue of a one-dimensional classical mechanical equation of motion in a 
time-dependent external potential. Moreover, as an additional approximation the actual 
substrate potential V{x, z) with a smooth lateral variation of Vatt (Eqs. (lA^) and (|A3|) ) can 



be replaced by the steplike potential Voo{x,z) = Q{—x)V-{z) + 6(a;)V"+(z). As it turns out 
both approximations have only negligible effects on the numerical results for the interface 
profiles (see Figs. ^ and ^ because for thick wetting films the width of the profile l{x) is 
much larger than the scale of the lateral variation of V{x, z). 

Figure |TT| shows liquid-vapor interface profiles for a SCS whose parameters are chosen 
such that first-order wetting transitions occur at T^ ^ 1.102 and T^ ~ 1.314 on w^ and 
i(;_, respectively. Here the liquidlike layers are very thin. Therefore also the width of the 
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profile is rather small as compared to a system undergoing a critical or complete wetting 
transition. In this case the replacement of the smooth substrate potential V{x, z) by the 
steplike potential Voo{x, z) yields a detectable difference but without changing the qualitative 
behavior. Although the applicability of the present sharp-kink approximation for such thin 
films certainly deserves a caveat we conclude that only in the case of a substrate potential, 
which gives rise to first-order wetting transitions, the transition region of the interface profile 
on structured substrates is confined to a region of the order of o"/ around the heterogeneity. 

B. Chemical step within a surface layer 

The previous subsection demonstrates that the properties of liquidlike wetting layers on 
a SCS heterogeneity are determined to a large extent by the asymptotic substrate potentials 
V±{z) of the semi- infinite substrate halves a; ^ 0. This is also the case for the LCS. The 
wetting properties of the substrate halves x ^ give rise to the asymptotic thickness of 
the wetting layer l{x — >• ±00) = l± as discussed in Subsec. [II C| and l{x) smoothly and 



monotonously interpolates between them. Figure |r^ displays a typical example for the 
influence of the thickness of the surface layer n on the profiles l{x). This system exhibits 
a critical wetting transition at Tj"* = T^* = T* = 1.2; we choose T* = 1.1 and A/i* = 
1.5 ■ 10~^. Without an overlayer (i.e., the number of inhomogeneous surface layers n = 0) 
the system corresponds to the homogeneous and planar substrate wh- As n is increased 
there is a crossover of l{x) to the profile of a SCS. The typical profile widths are - as for the 
SCS - given by Eqs. ([4.15|) and (|4.16|) if one of the substrate halves exhibits a complete or 
a critical wetting transition, respectively. In the case of a first-order wetting transition, i.e., 
for thin liquidlike layers the deviation of l{x) from its asymptotes is confined to a region of 
the order of cr/ around x = 0. 

From the analytical results for the VDWT one expects that the profile l{x) for a LCS 
approaches its asymptotes l± faster than in the case of a SCS, i.e., via 6l{\x\ -^ cxo) ~ x~^ 
(see Eqs. ([4.1CI|) - ([4.12|) ) rather than x~^ (Eqs. (^]8|) and ([4.9|) ). However, since as observed 



above the VDWT become relevant only for very large values of |x| without affecting the 
main deviation of l{x) from l± around x = 0, this more rapid decay is analytically accessible 
but does not play an important role for the overall behavior of l{x). 

The divergence of the function l{x) for x ^ cxd in the case that the substrate half a; > 
is completely and the other substrate half x < is only partially wet is given by Eqs. (|4.19 ) 



and ( [4. 231 ) for T > T+ and T = T+, respectively. This behavior is shown in Figs. |TB| and 



p!4| . The asymptotic behavior for x > of the profiles for small A/i and small T^ — T, 
respectively, is in excellent agreement with the numerical solutions of Eqs. ( |4.18| ) and ( [4.22[ ). 



C. Chemical stripe 

First we analyze a CST for which the homogeneous substrate w corresponding to the 
embedding substrate undergoes a critical wetting transition at T* = 1.2 whereas the homo- 
geneous substrate corresponding to the inserted slab Wst exhibits a critical wetting transi- 
tion at T* = 1.0. We consider thermodynamic states along the complete wetting isotherm 
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(T* = 1.1, A/i) so that at coexistence the substrate as a whole is only partially wet. The 
complete wetting transition on Wst is suppressed by the finite lateral extension of Wgt- 

Figure |15| displays liquid- vapor interface profiles for different values of A/x. The width of 
the transition regions at a; = ±a/2 increases as A/z becomes smaller. If this width is small 
compared with the width a of the stripe region the profile l{x) attains the equilibrium film 
thickness /+ on the stripe area. Therefore, if a is large the profile for a CST is approximately 
composed of two interfacial profiles corresponding to two SCS located at a; = —a/2 and 
X = +a/2. As the stripe width a is decreased for fixed values of T and A/i the region where 
l{x) attains the value /+ decreases and ultimately vanishes as a becomes smaller than the 



width of the transition region. This behavior is shown in Fig. 16 



Figure [T^ presents interface profiles for a CST with the chemical species of the slab and 
the surrounding substrate exchanged. Therefore the substrate exhibits a critical wetting 
transition at T* = 1.0 and it is already completely wet at coexistence as the temperature 
is raised towards T* = 1.2. The latter is the transition temperature of the critical wetting 
transition on the homogeneous substrate corresponding to Wgt- If we choose T* = 1.1 and 
decrease Afi the whole substrate undergoes a complete wetting transition and the "dent" 
induced by the existence of the stripe-shaped heterogeneity is smeared out. 

The above analysis demonstrates that the liquid leaks out of the chemical stripe if on the 
embedding material the fiuid is close to a complete or critical wetting transition. Figure p!8| 
shows that this leaking is absent if on the embedding material the fiuid can form only a 
thin wetting film. This situation prevails if the outer material leads to a first-order wetting 
transition at T~ and the temperature is chosen such that T < T~ . Figure |18| corresponds 
to a CST for which the slab as a homogeneous substrate undergoes a first-order wetting 
transition, too, at T+ with T+ < T < T~ . The confinement of the liquid to the stripe 
is achieved both off and at two-phase coexistence. In view of practical applications of 
chemically structured surfaces in the context of microfluidics this tells for which materials 
and for which thermodynamic states a good performance without leakage can be expected. 

In the strict sense of thermal equilibrium so far only a few solid- fiuid systems are known to 



exhibit a true first-order wetting phase transition at coexistence (e.g., He^ on Cs [55,66|); up 



to now critical wetting has been observed only on fiuid substrates (e.g., alkanes on aqueous 



solutions of salt or glucose p37| , |68[| ). The generic case is that fiuids wet substrates completely 
above the triple point |]50| , |69| . Although our theoretical framework is based on equilibrium 
statistical mechanics, the present model calculations for the morphology of wetting films 
can be applied directly also to nonvolatile liquids as long as their interaction with solid 
substrates can be modelled by appropriate effective interface potential u!±{l). Alternatively, 
these potentials can be inferred empirically from suitable experiments (see, e.g., Ref. [[7r!|,[7T[]). 
In this sense our conclusions are valid also for nonvolatile liquids adsorbed on "hydrophilic" 
or "hydrophobic" parts of an inhomogeneous substrate, for which a wide range of wettability 
characteristics, beyond the strict equilibrium conditions, can be arranged. 

The line tension Tcsrici) associated with a chemical stripe can be written as 

Tcsria) = a{a^^g - a^_g) + 2tscs + Sr^a) (5.2) 

where cr^^g is the surface tension of the homogeneous "-I-" or "— " substrate in contact with 
the vapor, tscs is the line tension associated with a single SCS, and 6T{a) with, as it turns 
out, 6T(a -^ oo) ~ a~^ is the effective interaction between the two line structures a distance 



^5 
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a apart from each other. tcst{o) generates a lateral force per unit length 

/(a) = ^^ (5.3) 

acting on the stripe and leading to a compression or dilation which is balanced by the elastic 
forces of the substrate, /(a) is defined such that a positive sign corresponds to a dilation 
of the stripe. Our numerical results show that for large stripe widths and independent of 
the wetting characteristics of the materials under consideration /(a) is practically constant 
with /(a) = cr^_g — (Jw+g + C'(a~^), i.e., in Eq. ( |5.2| ) for a > 4cr/ the first contribution 
is the dominating one (see Fig. |l^). For solid substrates the force is expected to have 
practically no effect because the compressibility of the substrate is very small. However, for 
fluid substrates such as fluids covered with Langmuir-Blodgett films the force can have a 
significant effect. Therefore we propose to test this effective force between line structures by 
adsorbing liquidlike wetting films on a liquid substrate decorated with a Langmuir-Blodgett 
film which contains a stripe of different material. The adsorbed wetting films will cause 
a change in the stripe width depending on the compressibility of the Langmuir-Blodgett 
films. Moreover, such experimental arrangements would facilitate to probe our predictions 
concerning critical wetting transitions, as reported in Refs. |]B^ and [^, on chemically 
structured substrates. 



VI. GENERAL MODEL FOR CHEMICALLY HETEROGENEOUS SUBSTRATES 

It is possible to generalize the methods developed within the previous sections to sub- 
strates with arbitrary heterogeneities. The only required input is the substrate potential 
which in general is a function V^(r) of all coordinates r = {r\\,z) = {x,y,z) and approxi- 
mately follows from the summation of pair potentials. The substrate potential ^(r), together 
with the fiuid-fiuid interaction, gives rise to an interface profile l{r\\) which in general de- 
pends on both lateral coordinates ry = {x,y) (see Fig. |20|). We assume that the substrate is 
fiat and located in the half space w = {r E M.^\z < 0}. We take into account the excluded 
volume at the substrate surface by introducing a spatially varying (i^(r||). The insertion of 
the sharp- kink density profile 

p(r) = Qiz - rf^(r|,))[0(/(r||) - z)pi + e(^ - /(r||))p,] (6.1) 

into Eq. ( |2.1| ) and the decomposition of Q[p] into bulk and subdominant contributions yields 

n[p{r)] = Afib(pg;T,/i) + a([/(r||)];r,/i; [w], [V]). (6.2) 

Concerning the expression for Qs we consider the thermodynamic limit and omit artificial 
contributions which stem from the truncation of the system. A = LxLyLz is the volume 
filled with fluid particles and Qb is the free energy density of the bulk vapor phase (see 
Eq. (p.7|)). The subdominant contribution in Eq. (|6.2|) reads 



ns[l{ri\)] = / dV||{Aa/(r||) + a^z(r||;d^(rj|)) + Si3(r||,[/(r||)]) 
Ja 
+c^(r,|, /(rii); rf^(r||)) + C(r||, /(r,,); K(ri|)])} (6.3) 
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where 



oo cc 

a^i{v\\-d^{v\\)) = -^ j dzt{z)+pi J dzV{r)-d^{r\\)nP (6.4) 

du-(rii) 

can be interpreted as a local, spatially varying wall-liquid surface tension. Equation ( |6.4| ) is 
a generalization of Eq. ( p.l7|) . T,ig{r^\, /(ry)) is the surface free energy density containing the 
cost in free energy for deforming the liquid-vapor interface: 



oo 



Kr||)-«(r||') 



^r^^lril, [/(r,))]) = ^19-^-^ j^ d\ ' Jdz I dz'w{\v^^ - r,, l\z- z'\) (6.5) 





within the nonlocal and 



sSrHr||,[/(r||)]) = a,,y/l + (V|,/(r|,))2 (6.6) 

within the local theory with Vy = [d^, dy). The latter expression is the leading term of the 

gradient expansion of E|^ . In Eqs. ( |6^ and (|6.6| ) S^g does not depend on the absolute 
value of 1{tc\\) measured from the substrate surface but only on the relative differences ^(ry) — 
/(rj,) so that aig is independent of /(ry). Such additional dependences are brought about by 
replacing the sharp-kink density profile in Eq. ( |6.1| ) by a smooth one whose tails are cut off 
by the surface (see, e.g., Refs. [^ and [^). As a generalization of Eq. ( pBf ) 



CX) oo 

cu(ry,Z(ry);rf^(ry)) = App, j dzt{z) ~ Ap f dzV{r) (6.7) 

Kr||)-du,(r||) Z(r||) 

is the "local" effective interface potential for the effective interaction between the substrate 
surface and the liquid-vapor interface. Finally, 



oo 
^2 



C{riiJirn);[dUr\\)]) = -Y fdzt{z)+piAp f dz 



t{z) 



d„(r||)-/(r||) 



~'~"o" / '^^'"11 ' / d^ dvw{\rn — ry '|, \v\) 



rfii,{r||)-d„{r|| ') z 

oo oo 



—piAp rf^^ii ' / dz dvw{\r\\—r\\'\,\v\) (6.8) 

dit,(r|| ')-«{r||) 2 

takes into account the effects due to the lateral variation of (i^(ry). C vanishes for (itu(ry) = 
diu = const. 

If the substrate contains inhomogeneities with large linear extensions, i.e., the substrate 
is translationally invariant along one of the lateral directions, the subdominant contribution 

23 



Qs decomposes into "true" surface and line contributions as discussed above for the SCS, 
the LCS, and the CST. In the absence of such an additional translational symmetry one 
cannot identify a genuine line contribution. 

The equilibrium liquid-vapor interface profile l{r\\) minimizes Qg- Inserting /(r|j) into Qg 
renders an the equilibrium, laterally varying wall- vapor surface tension cr^g(r|j) such that 



where (see Eq. (|0|)) 



/ d\\ or^,(r|j) = min n,[/(rjj)] = QMr\\)] (6-9) 

Ja {'(r|i)} 



+a;(r||,r(r||);rf^(r||))+C(r||J(r|j);[rf^(r||)]). (6.10) 



Within the local theory and for dw{r^\) = d^ Eq. ( |6.10|) reduces to 



or^g(rii) = Afife/(r||) + a^i{r\\) + a^^^l + (V||/(r||))2 + cj(r||, /(r,,)) (6.11) 

with the corresponding ELE 

V^(r||) \ , 9^(r||,/) 



'l + (V/(r||))2y ^^ «=/>!,) 

= Afi,-App;t(r(r|,)-rfJ+Apr(r|j,r(r||)). (6.12) 

This type of equation is often used (e.g., in Ref. |^) to study the macroscopic properties 
of liquidlike wetting layers and droplets on structured substrates. This section provides a 
microscopic basis for the underlying concept of a local, spatially varying surface tension as 
it is used in several studies (e.g., in Refs. [PT|-PB|]). 



VII. SUMMARY 

We have obtained the following main results: 

1. Based on the description of wetting phenomena on homogeneous substrates within the 
framework of density functional theory (Sec. ||, Fig. p, we have derived the systematic 
decomposition of the grand canonical potential of a fluid in contact with a chemically 
heterogeneous substrate into bulk, surface and line contributions. The minimum of 
the latter yields the equilibrium morphology of the adsorbed liquidlike wetting films 
and the associated line tension within mean field theory. As paradigmatic cases we 
have studied a simple chemical step (SCS, see Fig. |^ for n — i> oo and Subsec. [IIIA|) , a 
chemical step within a surface layer supported by a homogeneous substrate (LCS, see 
Fig. H and Subsec. |IIIB|) , and a chemical stripe (CST, see Fig. |] and Subsec. P!IIC|) . 
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2. Across a SCS the profiles of the hquid- vapor interface of the hquidhke adsorbed wetting 
film morphology interpolate between their asymptotic values corresponding to the 
wetting film thicknesses on the two individual, homogeneous substrates forming the 
SCS (Figs. ^ and ||). The curvature of the profiles changes sign near the position of 
the step (Fig. ^). On each side of the step the lateral width A± of this transition 
region is governed by the corresponding lateral correlation length C,\\ of the height- 
height correlation function (see Figs. ^ and |TU|, and Eqs. ( [4.15| ) and ( [4.1(j| )). This 
width diverges according to power laws which depend on whether a complete or critical 
wetting transition (Fig. ^ is approached. Near these transitions the dependence of the 
profiles on the lateral coordinate x and on ,^|| exhibits scaling properties (Eqs. ( |4.17|) 



and ([4.21|) ) with singular scaling functions (Fig. ^ such that at the wetting transitions 
the profiles diverge algebraicly as function oix (see Eqs. ( [4.19| ), ( [4.23| ), and Figs. |l3|and 
|l^). The corresponding amplitudes as well as those for the van der Waals tails of the 
profiles (see Eqs. ( |4.(j| ) and (O)) have been determined analytically. The various power 



laws are modified by retardation (Subsec. [IV E| ). Near first-order wetting transitions 



the interface profiles vary on a molecular scale (Fig. |TT|) and are sensitive to the details 
of the laterally varying substrate potential. 

The morphology of the wetting film across a chemical step within a surface layer (LCS) 
is similar to that on a SCS. The confinement of the chemical heterogeneity to a thin 
surface layer leads to more rapidly decaying van der Waals tails (Eqs. ( [4.1C )-( 4.12| )) 



as compared with the SCS and modified amplitudes for the power-law divergences 
as function of x (see the discussion after Eq. ([4.24|) ). The scaling functions for the 
SCS and the LCS are the same. The main structural features of the wetting film 
across a LCS are already induced by very few heterogeneous surface monolayers on 
the substrate (Fig. ^. 

4. If the chemical heterogeneity is confined to a stripe (CST) complete wetting of the inner 
region is inhibited by an incompletely wetted outer region (Fig. |l^). The leakage of 
the liquid into the outer region is governed by ^y corresponding to the latter. Thus 
the tight confinement of the liquid to the stripe can be accomplished by choosing an 
embedding outer material on which only thin wetting films can form. This confinement 
can be achieved even if the inner material prefers macroscopicly thick wetting films. 
In the latter case and for wide stripes with width a the interface profile has the shape 
of a semi-ellipse (Eq. (|4.25| )) and the excess coverage supported by the stripe scales as 
q3/2 (Y[g_ |X^_ The perturbation of the interface profile near one step by the distant 
step has been determined analytically (Eq. ( [4.26|) ) as well as the corresponding van 
der Waals tails characterizing the decay into the outer region (Eqs. (|4.13| ) and ([4.14|) ). 
Only for sufficiently wide stripes the film thickness on the stripe can attain the value 
/+ for the corresponding homogeneous case (Fig. 16). If, on the other hand, the outer 



region undergoes a complete or critical wetting transition but the inner region does 
not favor it, the thickening wetting film spills over into the stripe region and drags its 
wetting behavior along leaving behind a dent (Fig. |T^. The depth of the dent behaves 
nonmonotonously as function of the undersaturation Afi. Ultimately, for A/x -^ the 
depth of the dent vanishes and its width diverges proportional to ^\\ so that the net 
depletion of the coverage caused by the stripe diverges as (A/i)"'^/^^ (Fig. p!7|) . 
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5. The dependence of the hne tension rcsrici) associated with a chemical stripe on the 
stripe width a leads to a lateral force per unit length /(a) = —dTcsT{o)/da (Eqs. ( ^.21) 
and ( |5.3| )) acting on the stripe and leading to a compression or dilation which is 
balanced by the elastic forces of the substrate. Except for small stripe widths a ~ 
0"/ the constant contribution /o = cr^_g — cfw+g is the dominating one. The excess 
contribution fex{,o) = f{,o) — fo{a) decays as /ex(a -^ oo) ~ a~^ (Fig. |19D. Whereas 
for solid substrates this force has practically no effect we expect that for liquidlike 
wetting films on a liquid substrate decorated with a Langmuir-Blodgett film which 
contains a stripe of different material the force can become important. In this case 
the adsorbed wetting films cause a detectable change in the stripe width depending 
on the compressibility of the Langmuir-Blodgett films. 

6. In Sec. ^we have presented a systematic microscopic derivation of the Euler-Lagrange 
equation (Eq. ( |6.12| )) for a liquidlike film adsorbed on an arbitrary chemically struc- 
tured substrate (Fig. ^Oj). This approach also provides a microscopic calculus for deter- 
mining local, laterally varying wall-vapor and wall-liquid surface tensions (Eqs. ( |S.4| ) 
and ( p.lOp ) required for macroscopic analyses. Lateral inhomogeneities within the re- 
pulsive part of the substrate potential modify these expressions with respect to what 
is expected intuitively (see Eqs. ( |6.7| ) and ( |6.10| )). 
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APPENDIX A: POTENTIALS OF HETEROGENEOUS SUBSTRATES 

Approximately we determine the substrate potential by a pairwise summation of the 
interaction between a single fluid particle and all substrate particles. To this end we assume 
that the molecules in the substrate are located at orthorhombic lattice sites with the lattice 
constants gx, gy, and gz in the x, y, and z direction, respectively, which are taken to be con- 
stant throughout the substrate. The leading contribution from this summation corresponds 
to a three-dimensional integration over the substrate volume. The discrete sum generates 
in addition subdominant contributions which are proportional to powers of the lattice spac- 
ings. The substrate potential can be written as a sum V{x, z) = Vau{x, z) + Vrep{x, z) of an 
attractive and a repulsive contribution. In the attractive contribution we take into account 
the two leading orders of this power series whereas the repulsive contribution is modeled 
by a steplike crossover between the repulsive parts of the potentials of the corresponding 
homogeneous, semi-infinite substrates. This assumption is justified because the repulsive 
interaction decays very rapidly and is only significant for z < l-5o"/ where 2; = denotes the 
position of the nuclei of the top substrate layer. 

The different chemical species are distinguished by "+" and "— " denoting the different 
potential coefficients. If a homogeneous substrate is covered by a heterogeneous surface 
layer of finite thickness we denote its constituent molecules by "if" . For the interparticle 
potential we adopt the Lennard- Jones form 
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Under these conditions one finds for a simple chemical step (SCS), i.e., for two adjacent 
quarter spaces w+ and w_ (see Fig. || with n = oo) 
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with r = vxM-^ and 



Kr(^,-)=e(-a:)^ + 0(x)5 
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(A3) 



The coefficients -Ug , -^4 , and Wg are defined by the potentials of the respective homogeneous, 



fiat, semi-infinite substrates W7+ and W- (see Eq. ( p.l2|) ), whereas 



u 



4,x 



u^gxlgz- 



Similarly, for the homogeneous, fiat substrate wu covered by two different, adjacent 
surface layers ws^± which consist of n = d + 1 monolayers (compare Fig. 0) we obtain 
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where z = {z + dg^), r = yxM-^, r = \/x^+z^, and 
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LCS, 



V ^^(r z) — V 
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(A5) 



The coefficients u^ and u^ correspond to the coefficients of the homogeneous, fiat sub- 
strate covered by a homogeneous surface layer (see Eqs. ( |2.10| ) and ( |2.11| )) whereas u^'^ = 

H,± I 

M4 9x1 gz- 

Finally, for the slab w^t immersed in a homogeneous substrate w (see Fig. ^ the sum- 
mation yields 
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where the abbreviations X and F denote X = x — a/2 and F = x + a/2, and 
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with the coefficients defined as above for the SCS. 
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At large separations between a fiuid particle and the substrate retardation effects set in 
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APPENDIX B: DECOMPOSITION OF THE GRAND CANONICAL DENSITY 

FUNCTIONAL 

1. Single chemical step 

The sharp-kink ansatz (Eq. (p.3|)) for Eq. ( p.l| ) leads to the decomposition of the grand 
canonical free energy given by Eqs. ( p.4|) , ( p.7|) , and ( p.5|) with 

^^(4) = l±^^b + ^«,±« + (Tig + C^(4), (Bl) 

where the effective interface potentials uj±{l) are given by (compare with Eq. (|2.19|) ) 



I oo oo 

uo±{l) = I^p I Pi I dzt{z) - dzV±{z) 



(B2) 



Here we have omitted artificial contributions generated by truncating the system; these con- 
tributions are discussed in Ref . 0] . The surface contribution is that of two half substrates 
w^ and W- covered by liquidlike wetting layers of thickness /+ and /_ which are exposed 
to the potentials V^{z) and V^{z), respectively, of the corresponding semi-infinite, homo- 
geneous substrates. This defines a reference system such that the deviation of the smooth 
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profile l{x) from tlie reference configuration (see tlie daslied lines in Figs. ^ and |^) due to 
the heterogeneity and due to the smooth crossover of V{x,z) from V^{z) to V-^{z) leads to 
the line contribution 



ni[l{x)]=r{dt,h)+m^)]- 



(B3) 



The term T{d^,l±), which is independent of l{x), is given in Eqs. (B12), (B19), and (B20) 
in Ref. ||51|. (We note that there the + sign of the last term in Eq. (B20) must be replaced 
by a — sign.) The contribution which depends functionally on l{x) reads 

oo 

Oj[1{x)] = AD^bTex + / dx{uj{x, l{x); t/^) - uj{x, locix); rf+)} 
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oo oo l{x)-l(x') 
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w{x — X ,z — z) 



(B4) 



-oo — oo 



where loo{x) = /_G(— a;) + /+0(x). The first term measures the cost in free energy for 
replacing a certain volume of vapor by the liquid phase. It is proportional to the excess 
coverage 



dx{l{x) — loo{x)) 



(B5) 



and vanishes at coexistence /i = /iq. The second term is the integrated "local" effective 
interface potential 



oo oo 

uo{x, I] dw) = Ap ■^ p; / dz t{z) — dz V{x, z) ^ . 
The third term involving integrals of 

oo oo 

i{x,z) = / dx' / dz' w{x', z') 



(B6) 



(B7) 



takes into account the difference between d^ and d^ and vanishes for d^ = dtu- The last 
contribution with 



w[ 



{x,z)= I dyw{^/x^^+y^^+z^) 



(B8) 
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describes the free energy of the deformed free hquid-vapor interface. It is a nonlocal func- 
tional of l{x). A gradient expansion to the first order of this contribution yields the fourth 
term of the corresponding local functional uJioc[l{x)]: 



oo oo 



oo l{x)-l(x') 



-\(Apf j d. j dJ j d. I d/u,(|.-.1,|.-.'| 




(B9) 



2. Chemical stripe 

The sharp-kink approximation of the density profile for the CST is 

p{x, z) = {e (|x| - ^) e{z - rf-) + e (I - |x|) e{z - d^)} 

X [0(/(x) - z)pi + Q{z - lix))pg] . (BIO) 

Insertion of p{x, z) into Eq. ( p?T| ) leads to the decomposition 

n{[p{x,z)]-T,p-[^,[V]) = Kn,{pg,T,p) + An,{l_-T,p-[wl[V]) 

+Lyni{[l{x)]-T,p-[wi[V]). (Bll) 

with Qft given by Eq. ( p.7|) . The surface contribution stems from the reference system which 
in the case of the CST is the structured substrate covered by a liquidlike layer of thickness 
/_ that is exposed to the potential V-{z): 

0,(/_) = fi;(/_) (B12) 

with VL-{1-) defined by Eqs. (0) and (i|). 
The line contribution for the CST is given by 

^i[l{x)] = T{dtl^) + u[l{x)] (B13) 

and represents the free energy associated with the deviation of l{x) from the asymptotic 
value /_, i.e., of the reference configuration. In Eq. (P13|) the first term does not depend on 

l{x): 
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with SV{x,z) = V{x,z) - V-{z) and d^{x) = Q{\x\ - a/2)d- + Q{a/2 - |x|)rf+. The 
interpretation of the different terms in Eq. ( |B14|) is analogous to the SCS (compare Ref. WW )- 
The second term in Eq. ( |B13| ) depends on l{x): 



cj[/(x)] = Ai^foFg^ + / dx{uj{x,l{x)]d^) — uj{x,l-]d^)] 



— oo 

oo a/2 i{x)-d:;;, oo 



-oo -a/2 l(x)~d+ 2 

oo oo oo l{x)-l(x') 

— -(Ap)^ dx dx' I dz I dz' w{x — x', z — z') (B15) 



with 



-oo — oo 



V^:' = / dx{lix)-l-). (B16) 



The interpretation of the first two terms is the same as for the SCS; the third term is 
generated by the difference of the excluded volumes at the surfaces of the inner and the 
outer region and vanishes if d^ = d~. The last term is the same as in Eq. (p^) . 

The corresponding local functional expression ojioc[l{x)] follows by replacing the last term 
in Eq. ( |B15| ) by the local functional given in Eq. (p9|). We note that for all three substrate 
types (SCS, LCS, and CST) the functional ^^[/(a;)] (and therefore also oJioc[l{x)]) exhibit the 
same structure if d'^ = d~. 



31 



REFERENCES 

[1] W. Hansen, J. P. Kotthaus, and U. Merkt, in Semiconductors and Semimetals, edited by 
M. Reed (Academic, London, 1992), Vol. 35, p. 279; P. Bonsch, D. Wiillner, T. Schrimpf, 
A. Schlachetzki, and R. Lacmann, J. Electrochem. Soc. 145, 1273 (1998); J. Wang, 

D. A. Thompson, and J. G. Simmons, J. Electrochem. Soc. 145, 2931 (1998). 
[2] D. W. L. Tolfree, Rep. Prog. Phys. 61, 313 (1998). 

[3] Y. Xia and G.M. Whitesides, Annu. Rev. Mater. Sci. 28, 153 (1998). 

[4] F. Burmeister, C Schafle, B. Keilhofer, C Bechinger, J. Boneberg, and P. Leiderer, 

Adv. Mater. 10, 495 (1998). 
[5] J. A. Rogers, Z. Bao, and L. Dhar, App. Phys. Lett. 73, 294 (1998). 
[6] J. Aizenberg, A. J. Black, and G. M. Whitesides, Nature 394, 868 (1998). 
[7] S. Dietrich, Fluids in contact with structured substrates, proceedings of the NATO-ASI 

on New approaches to old and new problems in liquid state theory: inhomogeneities and 

phase separation in simple, complex, and quantum fluids, held at Patti Marina, Italy, 

July 7-17, 1998, edited by C. Caccamo, J. P. Hansen, and G. Stell (Kluwer, Dordrecht, 

1999), p. 147. 
[8] A. Paterson and M. Fermigier, Phys. Fluids 9, 2210 (1997). 
[9] G. Wiegand, T. Jaworek, G. Wegner, and E. Sackmann, J. Colloid Interf. Sci. 196, 299 

(1997). 
[10] P. Zeppenfeld, V. Diercks, C. Tolkes, R. David, and M. A. Krzyzowski, App. Surf. Sci. 

130-132, 484 (1998). 
[11] A. Karim, J. F. Douglas, B. P. Lee, S. C. Glotzer, J. A. Rogers, R. J. Jackman, 

E. J. Amis, and G. M. Whitesides, Phys. Rev. E 57, R6273 (1998). 
[12] T. Pompe, A. Eery, and S. Herminghaus, Langmuir 14, 2585 (1998). 

[13] L. Rockford, Y. Liu, P. Mansky, T. P. Russell, M. Yoon, and S. G. J. Mochrie, Phys. 

Rev. Lett. 82, 2602 (1999). 
[14] H. Gau, S. Herminghaus, P. Lenz, and R. Lipowsky, Science 283, 46 (1999). 
[15] B. S. Gallardo, V. K. Gupta, F. D. Eagerton, L. I. Jong, V. S. Craig, R. R. Shah, and 

N. L. Abbott, Science 283, 57 (1999). 
[16] C. S. Chen, M. Mrksich, S. Huang, G. M. Whitesides, and D. E. Ingber, Science 276, 

1425 (1997). 
[17] M. Matsuzawa, K. Kobayashi, K. Sugioka, and W. Knoll, J. Colloid Interf. Sci. 202, 

213 (1998). 
[18] E. Delamarche, A. Bernard, H. Schmidt, B. Michel, and H. A. Biebuyck, Science 276, 

779 (1997). 
[19] H. Shi, W.-B. Tsai, M. D. Garrison, S. Ferrari, and B. D. Ratner, Nature 398, 593 

(1999). 
[20] J. B. Knight, A. Vishwanath, J. P. Brody, and R. H. Austin, Phys. Rev. Lett. 80, 3863 

(1998). 
[21] M. Grunze, Science 283, 41 (1999). 
[22] R. F. Service, Science 282, 399 (1998). 

[23] J. L. Wilbur, A. Kumar, E. Kim, and G. M. Whitesides, Adv. Mater. 6, 600 (1994). 
[24] A. Kumar, H. A. Biebuyck, and G. M. Whitesides, Langmuir 10, 1498 (1994). 
[25] R. J. Jackman, J.L. Wilbur, and G. M. Whitesides, Science 269, 664 (1995). 
[26] P. C. Hidber, W. Helbig, E. Kim, and G. M. Whitesides, Langmuir 12, 1375 (1996). 

32 



[27] J. R. Henderson, J. Phys.: Condens. Matter 11, 629 (1999). 

[28] Y. Pomeau and J. Vannimenus, J. Colloid Interf. Sci. 104, 477 (1985). 

[29] M. O. Robbins, D. Andelman, and J.-F. Joanny, Phys. Rev. A 43, 4344 (1991). 

[30] C. J. Boulter, Phys. Rev. E 57, 2062 (1998). 

[31] P. Lenz and R. Lipowsky, Phys. Rev. Lett. 80, 1920 (1998). 

[32] P. S. Swain and R. Lipowsky, Langmuir 14, 6772 (1998). 

[33] R. Lipowsky, P. Lenz, and P. S. Swain, MPI fiir KoUoid- und Grenzfiachenforschung 

Teltow preprint (1999). 
[34] C Bauer, S. Dietrich, and A. O. Parry, preprint cond-mat/ 9904063] (1999). 



[35] S. Nechaev and Y.-C. Zhang, Phys. Rev. Lett. 74, 1815 (1995). 

[36] D. Urban, K. Topolski, and J. De Coninck, Phys. Rev. Lett. 76, 4388 (1996). 

[37] T. W. Burkhardt, J. Phys. A: Math. Gen. 31, L549 (1998). 

[38] M. W. Cole and E. Vittoratos, J. Low Temp. Phys. 22, 223 (1976). 

[39] L. Lajtar and S. Sokolowski, J. Chem. Soc. Faraday Trans. 88, 2545 (1992). 

[40] C. Chmiel, K. Karykowski, A. Patrykiejew, W. Rzysko, and S. Sokolowski, Mol. Phys. 

81, 694 (1994). 
[41] P. Rocken, A. Somoza, P. Tarazona, and G. Findenegg, J. Chem. Phys. 108, 8689 

(1998). 
[42] L. J. Douglas Frink and A. G. Salinger, J. Chem. Phys. 110, 5969 (1999). 
[43] S. K. Nath, P. F. Nealey, and J. J. de Pablo, J. Chem. Phys. 110, 7483 (1999). 
[44] W. Gac, A. Patrykiejew, and S. Sokolowski, Surf. Sci. 306, 434 (1994); Thin Sohd Films 

298, 22 (1997). 
[45] M. Schoen and D. J. Diestler, Chem. Phys. Lett. 270, 339 (1997); Phys. Rev. E 56, 

4427 (1997). 
[46] H. Bock and M. Schoen, Phys. Rev. E 59, 4122 (1999). 

[47] S. Curtarolo, M. W. Cole, M. J. Bojan, and W. A. Steele, Phys. Rev. E 59, 4402 (1999). 
[48] J. Israelachvili, Intermolecular and surface forces, 2nd ed. (Academic, London, 1991). 
[49] S. Dietrich and M. Schick, Phys. Rev. B 33, 4952 (1986). 
[50] S. Dietrich, in Phase Transitions and Critical Phenomena, edited by C. Domb and 

J. L. Lebowitz (Academic, London, 1988), Vol. 12, p. 1. 
[51] W. Koch, S. Dietrich, and M. Napiorkowski, Phys. Rev. E 51, 3300 (1995). 
[52] R. Evans, Adv. Phys. 28, 143 (1979). 
[53] J. D. Weeks, D. Chandler, and H. C. Andersen, J. Chem. Phys. 54, 5237 (1971); 

H. C. Andersen, J. D. Weeks, and D. Chandler, Phys. Rev. A 4, 1597 (1971). 
[54] N. F. Carnahan and K. E. Starling, J. Chem. Phys. 51, 635 (1969). 
[55] P. G. de Gennes, Rev. Mod. Phys. 57, 827 (1985). 
[56] M. Schick, in Liquids at interfaces. Proceedings of the Les Houches Summer School 

Lectures, Session XLVIII, edited by J. Chavrolin, J. F. Joanny, and J. Zinn- Justin 

(Elsevier, Amsterdam, 1990), p. 415. 
[57] D. E. Sullivan and M. M. Telo da Gama, in Fluid interf acial phenomena, edited by 

C. A. Croxton (Wiley, New York, 1986), p. 45. 
[58] R. Evans, in Fundamentals of inhomogeneous fluids, edited by D. Henderson (Marcel 

Dekker, New York, 1992), p. 85; see Sec. V therein. 
[59] (a) M. Napiorkowski and S. Dietrich, Phys. Rev. B 34, 6469 (1986); (b) Europhys. Lett. 

9, 361 (1989); (c) S. Dietrich and M. Napiorkowski, Phys. Rev. A 43, 1861 (1991). 



33 



[60] S. Dietrich, in Phase Transitions in Surface Films 2, proceedings of the NATO-ASI (Se- 
ries B) held in Erice, Italy, 19-30 June 1990, edited by H. Taub, G. Torzo, H. J. Lauter, 

and S. C. Fain, Vol. B 267, p. 391 (1991). 
[61] M. Kagan, W. V. Pinczewski, and P. E. Oren, J. Colloid Interf. Sci. 170, 426 (1995). 
[62] C Bauer and S. Dietrich, Eur. Phys. J. B, in press (1999) (preprint |cond-mat/9812266D . 
[63] R. Lipowsky, Phys. Rev. B 32, 1731 (1985). 

[64] M. P. Nightingale, W. F. Saam, and M. Schick, Phys. Rev. B 30, 3830 (1984). 
[65] E. Cheng, M. W. Cole, J. Dupont-Roc, W. F. Saam, and J. Treiner, Rev. Mod. Phys. 

65, 557 (1993). 
[66] R. Hallock, J. Low Temp. Phys. 101, 31 (1995); and references therein. 
[67] K. Ragil, J. Meunier, D. Broseta, J. O. Indekeu, and D. Bonn, Phys. Rev. Lett. 77, 

1532 (1996); N. Shahidzadeh, D. Bonn, K. Ragil, D. Broseta, and J. Meunier, Phys. 

Rev. Lett. 80, 3992 (1998). 
[68] T. Pfohl and H. Riegler, Phys. Rev. Lett. 82, 783 (1999). 
[69] D. Beysens, in Liquids at Interfaces, proceedings of the Les Houches Summer School 

Lectures, Session XLVIII, edited by J. Chavrolin, J.-F. Joanny, and J. Zinn- Justin 

(Elsevier, Amsterdam, 1990), p. 499. 
[70] H. I. Kim, C. M. Mate, K. A. Hannibal, and S. S. Perry, Phys. Rev. Lett. 82, 3496 

(1999). 
[71] S. Herminghaus, K. Jacobs, K. Mecke, J. Bischof, A. Fery, M. Ibn-Elhaj, and 

S. Schlagowski, Science 282, 916 (1998). 
[72] M. E. Fisher and A. J. Jin, Phys. Rev. B 44, 1430 (1991); Phys. Rev. Lett. 69, 792 

(1992). 
[73] A. O. Parry and R. Evans, Mol. Phys. 78, 1527 (1993). 



34 



FIGURES 




z=d 



w 



FIG. 1. Sharp-kink approximation for the particle density distribution. At the position z = Iq 
of the Hquid-vapor interface there is a stephke variation of the particle density between the bulk 
liquid (pi) and the bulk vapor density (pg). At z = L^ the density is truncated in order to facilitate 
the thermodynamic limit. Due to the repulsion between the fluid and the substrate particles the 
particle density vanishes for z < dw giving rise to an excluded volume. 
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FIG. 2. Morphology l{x) of a liquidlike film which covers a flat, heterogeneous substrate (layer 
chemical step (LCS)). The substrate surface is located at z = 0. In general the substrate consists of 
a homogeneous part (denoted as H) which is covered by two adjacent, semi-infinite surface layers 
composed of different chemical species ("+" and "— ") which form a sharp chemical boundary at 
X = 0. The surface layer consists oi n = d+1 molecular monolayers with a lattice spacing g^ in the 
z direction. The limit n ^ oo corresponds to a substrate which exhibits a single "simple chemical 
step" (SCS). /=p = /(x ^ Too) are the equilibrium film thicknesses of the corresponding substrates 
composed of "ff" plus a laterally homogeneous surface layer ("+" or "— "). We assume that 
"ff" covered by a "+" layer is the stronger substrate favoring the adsorption of thicker liquidlike 
films, i.e., /+ > L. The two substrate halves are characterized by different excluded volumes d^ 
(dash-dotted lines). The system is translationally invariant in the y direction and truncated at 
z = Lz and x = ^Lx/2 in order to facilitate the proper thermodynamic limit. The long-dashed 
line corresponds to the sharp-kink approximation for the lateral profile. 
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substrate 
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FIG. 3. Morphology of a liquidlike film which covers a flat substrate exhibiting a "chemical 
stripe" (CST). A slab consisting of chemical species denoted by "+" is immersed in a homogeneous 
"— " substrate. Thus in top view a chemical stripe with two sharp chemical steps is formed. The 
slab extends from x = —a/2 to x = a/2, i.e., the stripe width is a. L = /(|x| — > oo) is the 
equilibrium thickness of the liquidlike wetting film corresponding to the homogeneous substrate 
composed of "— " particles. In contrast to the SCS and the LCS (compare Fig. ^), here both 
cases /+ > I- and /-|- < /_ - where 1+ corresponds to the film thickness on the homogeneous "+" 
substrate - have to be considered. The slab and the embedding substrate exhibit different excluded 
volumes d^. The system is translationally invariant in the y direction and truncated at z = L^ 
and X = ^Lx/2 in order to facilitate the proper thermodynamic limit. 
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FIG. 4. Typical example for the shape of a liquid- vapor interface of a liquidlike film covering 
a SCS or LCS (full line) and the corresponding contour line U{x, z) = (dashed line) which 
separates regions with positive interface curvature {K{x,z) = U{x,z)/aig > 0) from regions with 
negative curvature {K{x,z) < 0) (see Eqs. ( [4.1D and ([4.2D ). The intersection between both lines 
is the turning point of the interface profile. The interface profile shown here is calculated for 
T* = kBT/ef = 1.3 and A/j* = A/j/ef = 10~^ and is also displayed in Fig. 0. The substrate 
is a SCS with the interaction parameters chosen such that the system exhibits critical wetting 
transitions at T^* = 1.0 and T^* = 1.2 on w-^- and W-, respectively (compare Fig. |7|). Therefore 
along the complete wetting isotherm (T* = 1.3, A^u) the functions Ai^f,l + uj±{l) corresponding to 
the homogeneous and flat substrates w^ and w- each exhibit only one extremal point and there 
is only a single contour line U{x, z) = interpolating between these two extremal points. 
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FIG. 5. Incomplete wetting of the substrate half x < and complete wetting of the substrate 
half X > can be obtained by approaching a thermodynamic state (T+ < Tq < T^ , Afi = 0) on dif- 
ferent thermodynamic paths in the (T, n) plane. This is illustrated by the corresponding schematic 
behavior of the equilibrium film thickness l± . For reasons of simplicity the coexistence curve hq (T) 
is straightened out. (a) Upon approaching ho{Tq) along the isotherm with T+ < Tq < T^ 1+ 
diverges whereas I- remains finite, (b) If the temperature is increased along coexistence fi = fiQ 
as T approaches T^ /_|_ diverges and is macroscopicly large for T^ < T < Tq. /_ also increases 
but remains finite since T < Tq < T^ . Both thermodynamic paths lead to the same final state 
{Tq, ^q{Tq)) of the system and thus to the same interfacial profile 1{x;Tq, hq(Tq)). This is also true 
if the wetting transition at T^ is first order; in this case in (b) l^ would jump to infinity at T = T^ . 
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FIG. 6. Scaling functions / and g as defined by Eqs. ( 4.17 ) and ( 4.21 ) calculated for 
Lennard- Jones potentials (C = 6). The scaling function / governs the evolution of the interface pro- 
file if the thermodynamic state (T+ < Tq < T~ , A^ = 0) at two-phase coexistence is approached 
along an isotherm (T = Tq, Afj.) ; for this limiting thermodynamic state the substrate half x > of 
a SCS or LCS is completely wet whereas the substrate half x < is only partially wet (see Fig. ^a). 
Analogously, the evolution of the interfacial profile upon approaching a critical wetting transition at 
(To = r+ < T^ , A^ = 0) along a path at coexistence as shown in Fig. ^3 is described by the scaling 
function g. In the limit ?/ ^ one has f{y ^ 0) ~ y^'^ and g{y ^ 0) '^ y^'^. f{y) and g{y) attain 
the value 1 for y ^ oo according to /(y ^ oo) = 1 - 6(^/3 - 1){V3 + l)'^e~^+^ e'y + 0{e-'^y) 
and g{y^oo) = l- Qe''^ e'^ -h 0{e-'^y). 
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FIG. 7. Shape of the hquid-vapor interface profile on a SCS (see Fig. |2|, with n = oo) along an 
isotherm (T* = 1.3, A^*), i.e., for complete wetting of both substrates w+ and w^. The parameters 
are chosen such that both substrates individually exhibit critical wetting transitions, Wj^ at T^ = 1.0 

and 



and w^ at T^ = 1.2: d^ 
0.277efO"f, whereas 



2.079e/-a^ 



u; 



12.475e/cj), 



~9x/gz 



uA 



'f^f 



u. 



= O.SOOUj*". Therefore for the temperature chosen here at two-phase 



coexistence both substrate halves are completely wet. The profiles shown here remain practically 
unchanged if the local ELE in Eq. ( ^.10 ) is approximated further by the square-gradient expression 



in Eq. (5.1) or if the substrate potential V{x,z) is replaced by the steplike substrate potential 

Voo{x, z) = V+{z)Q{x) + V-{z)Q{-x). 
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FIG. 8. Widths A± of the interface profiles (0, A) shown in Fig. |7| compared with the lateral 
height-height correlation lengths Ci determined by l± (full and dashed line, Eq. ( 4.15 )). A± is 



defined as that value |x| at which l{x) starts to deviate from its respective asymptote l± by 10% 
of |/+ — /_| and measures the lateral extension of the transition region of the interface profile. 
For A^ -^ both A± and ^j, diverge according to (A^)"^'^ as expected for complete wetting 
transitions. 
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FIG. 9. Shape of the hquid- vapor interface profile of a Hquidhke film across a SCS on a 
thermodynamic path at coexistence A/x = with the temperature T* increasing towards the 
wetting transition temperature T^ = 1.0 characteristic for both substrates w±. The parameters 
are chosen such that the wetting transitions are continuous: u- and u- ^ are equal to uf as in 
Fig. ^ but with nj = u'^^ = 14.035ej(Ti and d^ are also as in Fig. ^. Therefore 1+{T) / 1-{T) 
although T+ = T~ . Upon approaching the wetting transition temperature both /+ and /_ diverge 
as t~^ with tyj = {Tyj — T)/T^. As in Fig. M the results are practically insensitive against using the 
square-gradient approximation or using a steplike varying substrate potential. We note that here 
the scale of the x axis is about three times larger than in Fig. 0. 
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FIG. 10. Widths IS.± of the interface profiles (0, A) corresponding to the profiles shown in 
Fig. |9| as compared with the lateral height-height correlation lengths ^n determined by l± (full 
and dashed line, Eq. (4.16)). A-t is defined as that value \x\ at which l{x) starts to deviate from 



its respective asymptote l± by 10% of |/^ 
expected for critical wetting transitions. 



/_|. Both A-t and ^,, diverge as tw for t^^, — > as 



5/2 
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FIG. 11. Liquid-vapor interface profiles across a SCS with the parameters chosen such that 
both substrates individually exhibit first-order wetting transitions, the substrate W- at T^ « 1.314 
(with d~ = af, Ug = 2.513e/(Tp uj = nj^ = 3.770efai, and Ug 



at T* PS 1.102 (with d+ = l.05a 



0.335e/cr^) and the substrate w+ 



/' "3 



3.710efCJ 



/^/' 



u 



A,x 



5.566ej(Tp and Ug 



0.876e/o-p 



The system is at two-phase coexistence A// = 0. The full and dashed lines correspond to the full 
smooth substrate potential V{x,z) and its steplike approximation Voo{x^z), respectively. In both 
cases the local ELE (Eq. ( p.lO[ )) and its square-gradient approximation (Eq. (5T)) yield practically 
the same results. 
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FIG. 12. Liquid-vapor interface profiles on a LCS for different numbers n of monolayers in 
the inhomogeneous surface layer, n = corresponds to the homogeneous substrate wh, whereas 
n = oo leads to a SCS with the substrate halves w+ and w-. dw, uf , and uf^ are chosen 
equal to d", u~, and u~^ as in Fig. ^, respectively, such that both substrate halves exhibit a 
common critical wetting transition at T+* = T~* = T^ = 1.2. The interaction parameters of the 
surface layers are n^ = l.luf , u^ = 1.05uf^, and Qx = Qz- They do not affect the order of the 
wetting transition and its transition temperature but they lead to a difference between 1+{T) and 
l^[T). The temperature is fixed at T* = 1.1, i.e., both substrate halves are only partially wet at 
coexistence A/^ = 0. Moreover, we choose A//* = 1.5 • 10~^. As the number of monolayers in the 
surface layer n is increased the interface profile for a SCS composed of Wj^ and W- evolves out of 
the flat profile for the homogeneous substrate wh- Remarkably already a few adsorbed monolayers 
have a pronounced effect on the wetting film. 
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FIG. 13. Liquid-vapor interface profiles across a LCS with the substrate half x > completely 
wet and the substrate half x < only partially wet. The parameters d^, uf , uf, and Qx/Qz are 
chosen as in Fig. |l2|, u^ = 0.75uf , and we consider n = 10 inhomogeneous monolayers. With this 
choice of parameters the substrate half x > exhibits a critical wetting transition at T^* = 1.2 
whereas the substrate half x < undergoes a weakly first-order wetting transition at T^* ~ 1.306 
with the thickness /_ of the liquidlike film of the order of lOo"/ in the vicinity of T^*- The 
temperature is fixed at T* = 1.3 such that at coexistence A/i = the substrate half x < remains 
partially wet. At coexistence l{x) diverges according to l{x — > oo) = ^compX^'^ (see Eqs. ( [1.19 ) 



and ( 4.20| )). As in the case of the SCS, the width A_|_ of the interface profile for x > diverges 



according to (A/u)~^' ^ for A^ -^ 0; A_ remains finite even at A^ = 0, for which it is of the order 
of ^,7 ~ lO^aj. The inset provides a magnified view of the region around x = 0. For A/i -^ and 
X ^ oo the profiles exhibit scaling and are governed by the same scaling function f{y) as the SCS 
(Eq. (IOtD). 
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FIG. 14. Liquid-vapor interface profiles across a LCS along a thermodynamic path at coexis- 
tence A/u = approaching the critical wetting transition temperature r+* = 1.2 of the substrate 
half X > 0. Here we choose the same values of d,, 



n^, and Qx/Qz as in Figs. 12 and 13, u^ = uf , 



and u~ = 0.92u^ . The latter choice induces a weakly first-order wetting transition on the sub- 
strate half x < at T^ ~ 1.216 with a thickness of the liquidlike wetting film of the order of lOcrj 
slightly below the transition temperature. In the limit T y T^ 1+{T) diverges as (T^ — T)~^ 
whereas l-iT) remains finite. At T = T^ one has the power-law behavior l{x — > oo) = ^cHtx^ 
(see Eqs. ( |4.23| ) and ( [4.24D ). As in the case of the SCS, the width A+ of the interface for x > 
diverges as tw with t^ = {T^ — T)/T^ -^ 0, whereas A_ remains finite even at tw = 0, where 
it is of the order of ^,7 ~ IQ'^af. The inset provides a magnified view of the region around x = Q. 
For t^„ ^ and x — > oo the profiles exhibit scaling and are governed by the same scaling function 



g{y) as the SCS (Eq. (4.21)). 
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FIG. 15. Liquid-vapor interface profiles on a CST (full lines). The slab Wst is symmetric 
around x = and forms a chemical stripe at the substrate surface with a width a/af = 4000. The 
parameters for the substrate potential are chosen as in Fig. |^ i.e., the whole substrate undergoes a 
critical wetting transition at T^ = 1.2. The critical wetting transition at T^ = 1.0 corresponding to 
the slab Wst is suppressed due to the finite lateral extension of the slab. The profiles correspond to 
the temperature T* = 1.1 so that at coexistence A^ = the inhomogeneous substrate as a whole 
is only partially wet. Upon approaching coexistence /i — > /io the interface profile broadens. The 
dashed lines indicate the equilibrium film thicknesses corresponding to the flat and homogeneous 
"+" and "— " substrate; /+ = oo for A/x = 0. If the width A_|_ of the transition regions near 
|x| = a/2 is small compared to a (as it is the case for A//* = 10~^) one has l{x) Ri /+ in the middle 
of the stripe. 
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FIG. 16. Change of the morphology of the liquid-vapor interface profiles on a CST upon a 
variation of the stripe width a. The profiles correspond to a/af = 800, 2000, 4000, and 6000 from 
the inner to the outer one. The temperature T and the interaction potentials are chosen as in 
Fig. |l^ and A/x* = 10~^, i.e., the interface profile for a = AOOOaj is the same as in Fig. |l5|. If the 
stripe width is decreased the region where l{x) attains the equilibrium film thickness /+ (indicated 
by the dashed line) shrinks and ultimately disappears. Only for a >> cq, where oq ~ 2A_|_(/_|_) is a 
characteristic crossover width, l{x = 0) attains the value /+. One has A+ ~ .^n'" (compare Fig. @) 
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in the limit A/x 
750a f. 



0, i.e., /_|_ — > oo (see Eq. ( 4.15| )). For the present system 
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FIG. 17. Morphology of liquid-vapor interface profiles on a SST with the same parameters 



as in Fig. 15 but with the chemical species of the stripe and the surrounding region exchanged 



such that the slab Wgt favors thinner liquidlike films than the embedding substrate. Since the 
temperature T* = 1.1 is above the wetting transition temperature corresponding to the outer 
region the substrate as a whole is completely wet at coexistence A/i = 0. Upon approaching 
coexistence the equilibrium film thickness /_ corresponding to the surrounding substrate diverges 
and the "dent" induced by the presence of the stripe is deepening but finally it is smeared out and 
vanishes. Upon decreasing A^ at first the depth /_ — l{x = 0) of the "dent" increases, reaches 
a maximum value {I- — l{0))max ~ 25(T/ at A/u* ~ 3 • 10~^, and finally vanishes as (A/i)^'^ for 



A/i \ 0. Since the width of the dent is governed by the lateral correlation length ,^| 
which diverges as (A^u)^^'^ the depletion of the coverage jPea 
stripe diverges as (A/i)"^'-'^^. 
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FIG. 18. Liquid-vapor interface profiles on a CST which undergoes a first-order wetting 
transition for a choice of temperature such that the substrate corresponding to Wst is completely 
wet ( "hydrophilic" ) whereas the embedding substrate is only partially wet ("less hydrophilic" ) . The 



parameters are chosen as in Fig. 11. Therefore the outer region undergoes a first-order wetting 
transition at T^* ~ 1.314 whereas the homogeneous substrate corresponding to Wst exhibits a 
first-order wetting transition at T+* ~ 1.102. The temperature is taken as T* = 1.2 so that 
T+ < T < T~ . The width of the stripe is a = 400(7/ and the chemical potential is varied. Since 
the equilibrium film thickness /_ = /(|x| — > oo) is of the order of af the liquidlike "channel" is 
confined to the "hydrophilic" stripe region without leaking. The shape of the liquidlike channel for 
A^ = is in good agreement with the semi-elliptic shape given by Eq. ( 4.25 ). If the stripe width 
a is increased, l{x = 0) and l{x = 0) — I- increase as a^'^ and the excess coverage supported by 



the stripe diverges as T^x ~ a 
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FIG. 19. Excess contribution fex{o) = —d5T{a)/da to the force per unit length which acts 
on a chemical stripe. The system is the same as in Fig. |l^ but with T* = 1.0 and at two-phase 
coexistence A/i = 0. The total force is /(a) = /o + fexia-), i-e., it is the sum of fexio) and the 



constant contribution /o = cr.. 



w^g 



(^w+q which in the present example is /o 



0.l37ef/aj. In the 



hmit a ^ oo fex{o.) decays as a . Both /o and fex are positive, leading to a dilation of the stripe. 



55 



d„(r„) 




FIG. 20. Sketch of a fluid film on a general chemical inhomogeneity. The potential y(r|j,^) 
of the arbitrarily structured substrate (with rn = {x,y)) gives rise to a laterally varying exclusion 
length dyu{rii) and a liquid-vapor interface shape described by /(rn). The surface (r||,2; = /(rn)) 
separates the liquidlike wetting layer (below) from the bulk vapor phase (above). The system is 
truncated at x = ±La;/2, y = ±Ly/2, and z = L^ in order to facilitate the proper thermody- 
namic limit. (This truncation gives rise to artificial surface and line tensions generated by these 
boundaries.) The chemical inhomogeneity of the substrate is not indicated. 
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